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Abstract 

The  majority  of  the  world  shoreline  is  experiencing  some  form  of  erosion, 
which  will  become  more  serious  if  the  mean  sea  level  rise  accelerates 
because  of  the  greenhouse  effect.  The  recent  increase  of  coastal  storm 
damage  necessitates  the  development  of  numerical  models  for  predicting 
the  damage  progression  and  breaching  of  beaches,  coastal  stone  structures 
and  earthen  levees  during  extreme  storms.  This  report  summarizes  the 
continuing  effort  to  improve  our  quantitative  understanding  of  beach 
morphology  and  structural  damage  progression  with  the  goal  to  develop 
simple  and  robust  models  that  are  suited  for  engineering  applications.  Our 
effort  for  the  last  to  years  has  produced  the  cross-shore  numerical  model 
CSHORE,  which  is  presently  limited  to  the  case  of  alongshore  uniformity. 
CSHORE  consists  of  the  following  components:  a  combined  wave  and 
current  model  based  on  time-averaged  continuity,  cross-shore  and 
longshore  momentum,  wave  action,  and  roller  energy  equations;  a  sediment 
transport  model  for  suspended  sand  and  bedload;  a  permeable  layer  model 
to  account  for  porous  flow  and  energy  dissipation;  empirical  formulas  for 
irregular  wave  runup,  overtopping  and  seepage;  and  a  probabilistic  model 
for  an  intermittently  wet  and  dry  zone  for  the  purpose  of  predicting  wave 
overwash  and  structural  damage  progression.  The  wet  and  dry  model, 
which  is  the  latest  addition  to  CSHORE,  is  calibrated  and  verified  using  our 
107  small-scale  tests  for  irregular  wave  overtopping  and  overflow  on  a  levee 
as  well  as  118  Dutch  tests  for  low-exceedance  wave  overtopping,  in  which 
velocities  and  water  depths  were  measured  on  the  crest  and  landward 
(inner)  slope  of  dikes.  Finally,  the  computer  program  CSHORE  is  explained 
to  facilitate  use  and  modifications,  if  necessary.  In  the  near  future,  CSHORE 
will  be  compared  with  wave  overwash  experiments  and  structural  damage 
progression  experiments. 
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1  Introduction 

A  sand  beach  with  a  wide  berm  and  a  high  dune  provides  storm  protection 
and  damage  reduction,  recreational  and  economical  benefits,  and  biological 
habitats  for  plants  and  animals.  Most  sandy  beaches  are  eroding  partly  due 
to  sea  level  rise.  Beach  nourishment  is  widely  adopted  to  maintain  a  wide 
beach  for  a  developed  coastal  community  if  a  suitable  beachfill  is  available 
in  the  vicinity  of  an  eroding  beach.  Empirical  methods  based  on  field  data 
have  been  developed  for  the  design  of  beach  fills  (Headquarters,  U.S.  Army 
Corps  of  Engineers  (HQUSACE)  2002).  The  design  of  the  cross-shore 
beachfill  profile  is  normally  based  on  the  concept  of  an  equilibrium  beach 
profile.  The  alongshore  spreading  of  the  beachfill  is  generally  predicted 
using  a  one-line  model  coupled  with  the  Coastal  Engineering  Research 
Center  (CERC)  formula  or  the  formula  by  Kamphuis  (1991)  for  the 
longshore  sediment  transport  rate.  These  simple  beachfill  design  methods 
have  been  criticized  and  a  number  of  more  process-based  models  have  been 
proposed.  However,  the  process-based  models  may  not  necessarily  be  more 
accurate  at  present. 

Sediment  transport  is  caused  by  the  combined  action  of  waves  and  currents. 
Our  capabilities  of  predicting  wave  and  current  fields  have  improved 
steadily  over  the  last  30  years.  However,  the  predictive  capability  of  sedi¬ 
ment  transport  on  beaches  has  not  improved  to  the  same  degree.  The  major 
reason  for  this  discrepancy  is  that  no  dynamic  equation  is  available  to 
describe  the  motion  of  a  large  number  of  sediment  particles.  Consequently, 
sediment  transport  models  are  essentially  empirical  and  dependent  on 
reliable  sediment  transport  data.  Unfortunately,  sediment  dynamics  on 
beaches  are  highly  complex  and  involve  wide  ranges  of  morphological  scales 
in  time  and  space.  Correspondingly,  available  sediment  transport  models 
have  become  more  complex  and  less  transparent.  We  have  tried  to  synthe¬ 
size  available  data  and  formulas  to  develop  simple  and  transparent  formulas 
for  the  cross-shore  and  longshore  transport  rates  of  suspended  sand  and 
bedload  on  beaches.  The  simple  formulas  need  to  include  basic  sediment 
dynamics  sufficiently  so  that  the  formulas  will  be  applicable  to  small-scale 
and  large-scale  laboratory  beaches  and  eventually  prototype  beaches. 
Furthermore,  the  morphological  model  should  be  computationally  efficient 
to  allow  calibration  and  verification  using  extensive  data  sets.  The 
hydrodynamic  input  required  for  the  morphological  model  should  be 
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limited  to  quantities  that  can  be  predicted  routinely  and  reliably.  These 
considerations  have  guided  our  development  of  the  cross-shore  numerical 
model  CSHORE  presented  in  this  report. 

Coastal  storm  damage  has  been  increasing  mostly  due  to  the  recent  growth 
of  coastal  population  and  assets  and  possibly  due  to  the  intensification  of 
hurricanes  caused  by  global  warming.  Coastal  structures  including  earthen 
levees  (dikes)  have  been  designed  conventionally  for  no  storm  surge 
overflow  and  minor  wave  overtopping  during  a  design  storm.  Empirical 
formulas  for  wave  overtopping  rates  are  used  for  a  preliminary  design 
where  EurOtop  Manual  (2007)  recommends  the  latest  formulas.  Physical 
model  testing  is  normally  conducted  in  a  wave  flume  or  basin  for  a  detailed 
design.  Various  numerical  models  have  also  been  developed  to  predict 
detailed  hydrodynamics  that  are  difficult  to  measure  even  in  a  laboratory 
(Kobayashi  and  Otta  1987;  Kobayashi  1999;  van  Gent  2001).  The  latest 
numerical  models  for  hydrodynamics  are  reviewed  by  Losada  et  al.  (2008) 
and  Neves  et  al.  (2008).  However,  our  improved  predictive  capabilities  for 
the  hydrodynamics  have  not  really  improved  our  predictive  capability  for 
damage  progression  partly  because  damage  to  a  coastal  structure  is 
cumulative  (Melby  and  Kobayashi  1998).  As  a  result,  a  performance  or  risk- 
based  design  of  a  coastal  structure  relies  on  empirical  formulas  for  damage 
(e.g.,  Kobayashi  et  al.  2003).  This  practical  difficulty  is  similar  to  that  for 
sediment  transport  on  beaches.  Alternatively,  the  computationally-efficient 
CSHORE  calibrated  with  extensive  data  sets  has  been  developed  for  the 
design  of  inclined  structures  with  relatively  small  wave  reflection  where 
damage  progression  models  for  stone  structures  and  earthen  levees  will  be 
developed  by  modifying  the  sediment  transport  model  in  this  report.  The 
eventual  goal  is  to  predict  the  performance  of  an  inclined  structure  located 
on  a  movable  bottom. 
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2  History  of  CSHORE  Development 

The  history  of  the  cross-shore  model  CSHORE  is  summarized  to  provide 
an  overview  of  CSHORE  and  acknowledge  a  number  of  graduate  students 
and  visiting  scientists  who  contributed  to  the  development  of  CSHORE. 
The  present  version  of  CSHORE  includes  the  various  capabilities  added  to 
the  initial  model  developed  in  1998.  The  different  stages  of  development 
are  summarized  in  the  following  where  the  detail  of  each  stage  can  be 
found  in  the  listed  publications. 

The  cross-shore  model  CSHORE  was  developed  initially  to  predict  the 
cross-shore  transformation  of  irregular  nonlinear  waves  using  the  time- 
averaged  continuity,  momentum,  and  wave  energy  equations  together  with 
a  non-Gaussian  probability  distribution  of  the  free  surface  elevation. 
However,  empirical  formulas  of  limited  generality  were  required  to 
parameterize  the  wave  nonlinearity.  The  present  version  of  CSHORE  is 
based  on  linear  wave  theory  and  the  Gaussian  probability  distribution  to 
reduce  the  degree  of  empiricism  (Kobayashi  et  al.  1998;  Kobayashi  and 
Johnson  1998;  Johnson  and  Kobayashi  1998;  Kearney  and  Kobayashi 
2000;  Johnson  and  Kobayashi  2000). 

The  next  stage  of  the  CSHORE  development  was  motivated  by  the  need  of 
a  computationally-efficient  time-averaged  model  that  can  be  used  for  the 
design  of  porous  coastal  structures.  The  linear-wave  version  of  the  initial 
CSHORE  was  modified  to  account  for  the  effects  of  a  permeable  layer  for 
the  case  of  normally  incident  waves.  The  permeable  version  of  CSHORE 
was  called  CSHOREP.  The  permeability  effects  are  extended  to  obliquely 
incident  waves  in  the  present  CSHORE  (Meigs  and  Kobayashi  2004; 

Meigs  et  al.  2004;  de  los  Santos  and  Kobayashi  2005;  Ota  et  al.  2006; 
de  los  Santos  et  al.  2006;  de  los  Santos  and  Kobayashi  2006;  Kobayashi 
et  al.  2007b;  Kobayashi  and  de  los  Santos  2007;  Ota  et  al.  2007; 

Kobayashi  et  al.  2008b). 

Concurrently,  the  impermeable  version  of  CSHORE  is  being  extended  to 
predict  the  cross-shore  and  longshore  transport  rates  of  suspended  sand 
and  bedload  on  beaches.  A  series  of  extensions  were  made  in  the  following 
publications  to  make  the  model  both  versatile  and  better  verified  (Zhao  and 
Kobayashi  2005;  Kobayashi  et  al.  2005;  Agarwal  and  Kobayashi  2005; 
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Schmied  et  al.  2006a  and  b;  Agarwal  et  al.  2006;  Payo  et  al.  2006; 
Kobayashi  et  al.  2007a;  Buck  et  al.  2007;  Gencarelli  et  al.  2007;  Kobayashi 
et  al.  2008a  and  b;  Payo  et  al.  2009  Buck  et  al.  2008;  Gencarelli  et  al. 

2008a  and  b.  The  impermeable  and  permeable  versions  of  CSHORE  have 
been  merged  in  the  present  CSHORE  to  expand  the  range  of  practical 
applications. 

The  following  papers  summarized  the  progress  of  the  CSHORE  develop¬ 
ment  concisely  up  to  20o8(Kobayashi  2006;  Kobayashi  2008;  Kobayashi 
et  al.  2008e).  These  publications  were  based  on  the  earlier  version  of 
CSHORE  limited  to  the  wet  zone  below  the  mean  water  level.  To  extend 
CSHORE  to  the  zone  which  is  intermittently  wet  and  dry,  laboratory 
experiments  were  conducted  for  wave  overtopping  and  overflow  on  fixed 
levees.  The  laboratory  data  were  used  for  the  development  of  a  probabilistic 
model  for  the  wet  and  dry  zone  on  an  impermeable  bottom.  This  hydro- 
dynamic  model  coupled  with  the  sediment  transport  model  in  CSHORE  has 
been  used  to  predict  wave  overwash  of  dunes.  The  hydrodynamic  model  has 
also  been  extended  to  the  wet  and  dry  zone  on  a  permeable  bottom  for  the 
prediction  of  wave  overtopping  of  rubble  mound  structures.  This  model 
coupled  with  the  CSHORE  bedload  formula  modified  for  stone  has  been 
shown  to  be  capable  of  predicting  the  evolution  of  damaged  stone  armor 
layers  (Farhadzadeh  et  al.  2007;  Kobayashi  et  al.  2007c;  Farhadzadeh  et  al. 
2008;  Kobayashi  and  Farhadzadeh  2009;  Johnson  et  al.  2009;  Figlus  et  al. 
2009;  Figlus  et  al.  2010;  Farhadzadeh  et  al.  2009;  Farhadzadeh  et  al.  2010; 
Kobayashi  et  al.  2010a;  Kobayashi  et  al.  2010b). 
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3  Wave  and  Current  Models 

Cross-shore  sediment  transport  on  beaches  has  been  investigated 
extensively  (e.g.,  Kriebel  and  Dean  1985;  van  Rijn  et  al.  2003)  but  we  still 
cannot  predict  beach  profile  evolution  accurately.  To  improve  our  predictive 
capabilities,  sediment  transport  models  have  become  more  sophisticated 
but  less  transparent.  For  example,  Thornton  et  al.  (1996)  and  Gallagher 
et  al.  (1998)  used  the  energetics-based  total  load  model  of  Bailard  (1981)  to 
explain  the  offshore  movement  of  a  bar  at  Duck,  NC,  during  storms.  The 
onshore  bar  migration  on  the  same  beach  was  predicted  by  both  Hoefel  and 
Elgar  (2003),  using  the  skewed  acceleration  effect  on  bedload,  and 
Henderson  et  al.  (2004),  using  a  suspended  sediment  model.  The  contribu¬ 
tion  of  bedload  and  suspended  load  are  not  clear  at  present.  Kobayashi  et  al. 
(2008a)  made  an  attempt  to  synthesize  and  simplify  existing  cross-shore 
sediment  transport  models  with  the  aim  of  developing  a  simple  and  robust 
model  that  is  suited  for  engineering  applications  including  the  berm  and 
dune  erosion.  This  model  has  been  extended  to  predict  the  cross-shore  and 
longshore  transport  rates  of  bedload  and  suspended  load  under  the 
combined  wave  and  current  action  predicted  by  the  time-averaged, 
probabilistic  model  by  Kobayashi  et  al.  (2007a). 

Sediment  transport  on  beaches  is  caused  by  the  combined  action  of  waves 
and  currents.  The  hydrodynamic  input  required  for  a  sediment  transport 
model  depends  on  whether  the  sediment  transport  model  is  time- 
dependent  (phase-resolving)  or  time-averaged  over  a  number  of  waves.  A 
time-dependent  sediment  transport  model  such  as  that  by  Kobayashi  and 
Johnson  (2001)  is  physically  appealing  because  it  predicts  intense  but 
intermittent  sand  suspension  under  irregular  breaking  waves  (Kobayashi 
and  Tega  2002).  However,  the  time-dependent  model  requires  considerable 
computation  time  and  is  not  necessarily  more  accurate  in  predicting  slow 
morphological  changes  than  the  corresponding  time-averaged  model 
presented  in  the  following.  Horizontally  two-dimensional  (2-D)  wave  and 
current  models  are  presented  first  before  the  cross-shore  model  CSHORE, 
which  is  based  on  the  assumption  of  alongshore  uniformity. 

Figure  1  shows  obliquely  incident  irregular  waves  on  an  essentially  straight 
shoreline  where  the  cross-shore  coordinate  x  is  positive  onshore  and  the 
longshore  coordinate  y  is  positive  in  the  down  wave  direction.  The  beach  is 
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Figure  1.  Definition  sketch  for  incident  irregular  waves  and  wind  on  beach. 


assumed  to  be  impermeable.  The  depth-averaged  cross-shore  and  longshore 
velocities  are  denoted  by  U  and  V,  respectively.  Incident  waves  are  assumed 
to  be  unidirectional  with  0  equal  to  the  incident  angle  relative  to  the  shore 
normal.  The  height  and  period  of  the  irregular  waves  are  represented  by  the 
root-mean-square  (rms)  wave  height  Hrms  and  the  spectral  peak  period  Tp 
specified  at  the  seaward  boundary  located  at  x  =  o.  The  wave  angle  0  is 
assumed  to  be  in  the  range  of  1 0 1  <90  degrees  to  ensure  that  the  incident 
waves  propagate  landward.  The  wind  speed  and  direction  at  the  elevation  of 
10  m  above  the  sea  surface  are  denoted  by  W10  and  Qw,  respectively. 

The  mean  water  depth  h ,  with  the  over  bar  indicating  time-averaging  ,is 
given  by: 


h  =  (77-j-S  — zb)  (1) 

where  rj  =  wave  setup  above  the  still  water  level  (SWL) ;  and  S  =  storm 
tide  above  the  datum  Zb  =  o  which  is  assumed  to  be  uniform  in  the 
computation  domain,  where  Zb  =  bottom  elevation  relative  to  the  datum 
z  =  o  with  z  =  vertical  coordinate  taken  to  be  positive  upward.  Linear  wave 
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and  current  theory  for  wave  refraction  (e.g.,  Phillips  1977;  Mei  1989; 
Dalrymple  1988)  is  used  to  predict  the  spatial  variations  of  Hrms  and  0.  The 
dispersion  relation  for  linear  waves  is  expressed  as: 

co2  =  kgtanh([khy,a)p  =o)  +  k(KQxcos6  +  QlJsin6^j/ h  (2) 

where  co  =  intrinsic  angular  frequency;  k  =  wave  number;  g  =  gravitational 
acceleration;  cop  =  absolute  angular  frequency  given  by  cop  =  2jt  /  Tp ;  Qx  and 

Qy  =  time-averaged  volume  flux  per  unit  width  in  the  x  and  y  directions, 
respectively,  and  0  =  incident  wave  angle.  Equation  2  can  be  solved  itera¬ 
tively  to  obtain  k  and  co  for  known  cop,  h,  0,  Qx  and  Qy .  The  phase  velocity 

C  and  the  group  velocity  Cg  are  given  by: 


C  =  a/k;Cg 


nC;n 


2  kh 

sinh^2khj 


(3) 


The  wave  angle  0  is  computed  using  the  irrotationality  of  the  wave  number: 

J^-(ksin0)-J^-(ksin0)  =  0  (4) 

The  rms  wave  height  Hrms  defined  as  Hrms  =  \[8on  with  on  =  standard 

deviation  of  the  free  surface  elevation  ( rj)  is  computed  using  the  wave 
action  equation: 


d 

E 

(  0  ) 
C  ms  04--^- 

d 
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E 

f  Q  ' 
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dx 
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CO 

with 


E  =  pGo;  =  |p  gHl,  (6) 

where  E  =  specific  wave  energy;  p  =  fluid  density;  and  Db  and  Df  =  wave 
energy  dissipation  rate  per  unit  horizontal  area  due  to  wave  breaking  and 
bottom  friction,  respectively.  The  formulas  for  Db  and  Df  are  presented 
later  in  relation  to  the  cross-shore  model  CSHORE. 


ERDC/CHL  TR-12-22 


8 


The  time-averaged  volume  fluxes  Qx  and  Qy  in  Equation  2  are  expressed  as: 

Q<  =  hU  +  Qm;Q:J  =  hV  +  Qm  (7) 


with 


oo2cos0 


(8) 


where  U  and  V  =  time-averaged,  depth-averaged  velocities  in  the  x  and  y 
directions;  Qwx  and  Qwy  =  wave-induced  volume  fluxes  in  the  x  and  y 
directions;  (ga2/C)  =  volume  flux  due  to  linear  waves  propagating  in  the 

direction  of  0;  and  qr  =  volume  flux  of  a  roller  on  the  front  of  a  breaking 
wave.  The  roller  volume  flux  qr  is  estimated  using  the  roller  energy 
equation,  as  explained  by  Kobayashi  et  al.  (2005,  2007a): 

J-(pC2qr  cos0)  +  J^(p  C2qr  sin0)  =  DB-  Dr  (9) 


with 


Dr  =P9p,.Qr ;  P, .  =(0.1  +  Sft)>0.1 

(10) 

f)z,  _  87, 

S,  =  — 2_cos  0H - 2_sin0 

(11) 

b  dx  dy 

where  Dr  =  roller  dissipation  rate;  j 6r  =  wave-front  slope;  Si,  =  bottom  slope 
in  the  direction  of  wave  propagation.  The  wave  front  slope  j 6r  is  assumed  to 
be  0.1  unless  it  is  increased  by  the  positive  bottom  slope  Sb. 

The  mean  water  depth  h  and  the  current  velocities  U  and  V  are 
computed  using  the  time-averaged  continuity  and  momentum  equations 
(Phillips  1977;  Svendsen  et  al.  2002): 


(12) 
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(16) 


Svv  =  (n£  +  M  )cos20  +  £ 


1 

n - 

2 


Mr=pCqr ; 


(17) 


STy  =  (n£  +  Mr)cos0sin0  ;  Syy  =  (nE  +  Mr)  sin2  0  +  E 


1 

n - 

2 


(18) 


where  nx  and  ny  =  bottom  shear  stresses  in  the  x  and  y  directions;  tsx  and 
Tsy  =  wind  stresses  on  the  sea  surface  in  the  x  and  y  directions;  and  SW,  Sxy 
and  Syy  =  radiation  stresses  including  the  momentum  flux  Mr  of  a  roller 

o 

propagating  with  the  phase  speed  C.  It  is  noted  that  the  terms  Q  ,Q  Q 

LUX  LUX  LU  LJ 
o 

and  Qwi/  in  Equations  15  and  16  included  by  Phillips  (1977)  are  of  4-th  order 

in  terms  of  the  wave  height  and  are  normally  neglected.  The  present  circula¬ 
tion  model  based  on  Equations  12-18  is  a  simplified  version  of  SHORECIRC 
(Svendsen  et  al.  2002)  for  irregular  waves  where  SHORECIRC  assumes 
monochromatic  waves.  The  formulas  for  nx,  ny,  rsx  and  Tsy  are  presented  later 
in  relation  to  the  cross-shore  model  CSHORE. 


A  horizontally  2-D  model  C2SHORE  is  also  being  developed  (Johnson  and 
Grzegorzewski,  2011).  The  directional  spectral  wave  model  STWAVE  (Smith 
et  al.  2001)  is  used  to  predict  the  wave  transformation.  The  wave-induced 
fluxes  Qwx  and  Qwy  and  the  radiation  stresses  Sxx,  Sxy  and  Syy  are  computed 
from  the  predicted  directional  wave  spectra.  The  roller  effects  included  in 
Equations  8, 17,  and  18  are  neglected.  The  circulation  model  is  based  on 
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Equations  12-16  with  the  formulas  for  ti,x,  t by,  tsx  and  rsy  used  in  CSHORE. 
The  wave  and  circulation  models  are  coupled.  The  wave  field  is  computed  to 
estimate  twx  and  rwy  given  by  Equations  15  and  16  for  the  circulation  model 
which  computes  the  wave  setup  and  wave-induced  currents.  An  efficient 
finite  difference  method  is  used  to  solve  Equations  12-14  as  described  in  Shi 
et  al.  2007.  Johnson  and  Grzegorzewski  (2011)  compared  C2SHORE  results 
with  morphologic  change  data  for  Ship  Island  in  Mississippi  Sound,  with 
good  results. 
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4  Combined  Wave  and  Current  Model  in  Wet 
Zone 

The  cross-shore  model  CSHORE  assumes  alongshore  uniformity  but 
computes  the  wave  and  current  fields  simultaneously.  The  depth-integrated 
continuity  equation  of  water  given  by  Equation  12  requires  that  the  cross¬ 
shore  volume  flux  Qx  is  constant  and  equal  to  the  wave  overtopping  rate  q0 
at  the  landward  end  of  the  computation  domain.  Equations  7  and  8  yield: 

2 

Qx=hU  +  ^—!LcosQ  +  qrcosQ  =  q0  (19) 

go2 

Qy  =  hR  +  — Isin0  +  qrsin0  (20) 

C 

where  h  =  mean  water  depth;  U  =  mean  cross-shore  velocity;  which  is 
negative  and  offshore  because  cos  0  >  o  if  q0  =  o  (no  wave  overtopping); 
g  =  gravitational  acceleration;  an  =  standard  deviation  of  the  free  surface 
elevation  r\;  C  =  linear  wave  phase  velocity  in  the  mean  water  depth  h 
corresponding  to  the  spectral  peak  period  Tp;  and  qr  =  volume  flux  of  a 
roller  on  the  front  of  a  breaking  wave.  If  the  incident  wave  angle  0  is  small, 
Equation  20  can  be  approximated  by  Q  -hV  for  most  applications. 

For  the  case  of  alongshore  uniformity,  Equation  4  reduces  to  Snell’s  law: 

k  sin  0  =  constant  (21) 

which  is  used  to  obtain  the  wave  direction  0. 

The  constant  value  is  obtained  from  the  values  of  0,  h ,  and  Tp  specified  at 
the  seaward  boundary  x  =  o  located  outside  the  surf  zone  where  co  can  be 
approximated  by  cop  in  Equation  2.  Reflected  waves  are  neglected  in  this 
model. 

The  cross-shore  and  longshore  momentum  Equations  13  and  14  are 
simplified  as: 
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(22) 


(23) 


where  Sxx  =  cross-shore  radiation  stress;  p  =  water  density;  t bx  =  cross¬ 
shore  bottom  stress;  tsx  =  cross-shore  wind  stress  on  the  sea  surface;  Sxy  = 
shear  component  of  the  radiation  stress;  ny  =  longshore  bottom  stress; 
and  Tsy  =  longshore  wind  stress  on  the  sea  surface.  The  wind  shear  stresses 
may  not  be  negligible  especially  outside  surf  zones  on  natural  beaches 
(Lentz  et  al.  1999).  Linear  wave  theory  for  progressive  waves  is  used  to 
estimate  Sxx  and  SXy  as  in  Equations  17  and  18.  For  the  sake  of  clarity,  the 
equations  are  again  given  as: 


=  (nE  +  M„)  cos2  Q  +  E 


1 

n - 

2 


;  S/  —  [nE  +  M,.)cos0sin0  (24) 


with 

n  =  Cg/C;E  =  p ga\ ;  Mr  =  p Cqr  (25) 


where  Cg  =  linear  wave  group  velocity;  E  =  specific  wave  energy  with  the 
rms  wave  height  defined  as  Hrms  =  V8  an;  and  Mv  =  momentum  flux  of  a 
roller  propagating  with  the  phase  velocity  C. 

The  time-averaged  bottom  shear  stresses  in  Equations  22  and  23  are 
written  as: 


=|p fbuu„  ;t„„  =ip/„W/0  ;£/„  ={u2+V2fS  (26) 

where  U  =  depth-averaged  cross-shore  velocity;  V  =  depth-averaged 
longshore  velocity;/^  =  bottom  friction  factor;  and  the  overbar  indicates 
time  averaging.  The  bottom  friction  factor  fi,  is  of  the  order  of  0.015  but 
should  be  calibrated  using  longshore  current  data  because  of  the  sensitivity 
of  longshore  currents  to/5.  The  equivalency  of  the  time  and  probabilistic 
averaging  is  assumed  to  express  nx  and  ny  in  terms  of  the  mean  and 
standard  deviation  of  the  depth-averaged  velocities  U  and  V  expressed  as: 
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U  =  o  F  -1 7  —  n  F  •  U  =  o  F  •  F  =  F  +F  • 

u  'JT1U’V  'JT1V’ua  'JT1a’1a  X1  U  '  1  V  )  » 


0.5 


(27) 


with 


U 


V 


Fu  =  U*  +r  cosG  ;Fv=V,  +rsin0 ;  C7.  =  — ;  K  =  — ; 


a7 


a7 


(28) 


where  U  and  V  =  depth-averaged  cross-shore  and  longshore  currents; 

<jt  =  standard  deviation  of  the  oscillatory  (assumed  Gaussian)  depth- 
averaged  velocity  Ut  with  zero  mean;  and  r  =  Gaussian  variable  defined  as 
r  =  Ut/gt  whose  probability  density  function  is  given  by: 


(  o  \ 


f(r) 


exp 


r 

2 


(29) 


Linear  progressive  wave  theory  is  used  locally  to  express  Ut  in  terms  of  the 
oscillatory  free  surface  elevation  (r|  —  r|) : 


(30) 

which  yields  the  standard  deviation  arof  the  oscillatory  velocity  Ut-. 

oT  =  Co* ;  o,  =  o^  /  h  (31) 

It  is  noted  that  that  U»  =U  /  a T  and  V,  =  V/oT  are  of  the  order  of  unity  or 
less.  The  standard  deviations  of  U  and  V  are  given  by: 

au  =oTcosQ  ;ov  =  or|sin0|  (32) 

where  cos  0  >  o  but  sin  0  can  be  negative.  Substitution  of  Equation  27  into 
Equation  26  yields: 


p 

=  -pfiUrGhx  ;  xby 


=  \pfbv2TGby 


(33) 


with 
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Gbx  =  f  FvFj{r)dr  ;Gby  =  f  FvFaf(r)dr  (34) 

— oo  — oo 

which  must  be  integrated  numerically. 

The  wind  shear  stress  in  Equations  22  and  23  are  expressed  as: 

=  PaCDKo  cos6w  ;  Tsy  =  p aCDW?0  sinG^  (35) 


where  pa  =  air  density  (pa  =  1.225  kg/m3);  Cd  =  drag  coefficient,  W10  = 
10-m  wind  speed;  and  Qw  =  wind  direction  defined  in  Figure  1.  The 
formula  by  Large  and  Pond  (1981)  is  used  to  estimate  Cd  where  Cd  = 
0.0012  for  W10  <  11  m/s  and  Cd  =  (0.00049  +  0.000065  Wio)  for  W10  > 
11  m/sec.  It  is  noted  that  the  measured  values  of  Cd  during  tropical 
cyclones  by  Powell  et  al.  (2003)  indicated  no  increase  of  Cd  with  the 
increase  of  W10  much  above  25  m/sec.  In  short,  available  data  are 
insufficient  to  estimate  Cd  for  extreme  wind  conditions. 


The  wave  action  Equation  5  for  the  case  of  alongshore  uniformity  becomes: 


d_  E_ 
dx  co 


Cffcos9  + 


DB  +  F>f 
co 


(36) 


which  reduces  to  the  wave  energy  equation  if  co  is  constant  and  Qx  =  o.  In 
this  case,  the  wave  energy  equation  is  given  by: 


Db 


Df  ;FX  —  ECg  cosG 


(37) 


where  Fx  =  cross-shore  energy  flux  based  on  linear  progressive  wave 
theory;  and  Db  and  Df  =  energy  dissipation  rates  due  to  wave  breaking  and 
bottom  friction,  respectively.  The  energy  dissipation  rate  Db,  due  to  wave 
breaking,  in  Equation  36  is  estimated  using  the  formula  by  Battjes  and 
Stive  (1985),  which  was  modified  by  Kobayashi  et  al.  (2005)  to  account  for 
the  local  bottom  slope  and  to  extend  the  computation  to  the  lower  swash 
zone.  The  modified  formula  is  expressed  as: 


db  = 


= 


p  gaQHl  O- 1 


4  T 
0.88 


tanh 


’HQ) 

y  kh  ' 


H„ 


H„ 


0.88 


;a 


2  nS± 
3kh 


>1 


(38) 
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where  a  =  slope  effect  parameter;  Q  =  fraction  of  breaking  waves;  Hb  = 
breaker  height  used  to  estimate  Db;  T  =  intrinsic  wave  period  given  by 
T  =  27t/co  with  co  obtained  using  Equation  2;  Hrms  =  V8on  =  local  rms  wave 

height;  Hm  =  local  depth-limited  wave  height;  k  =  wave  number;  h  =  mean 
water  depth  including  wave  setup;  y  =  empirical  breaker  ratio  parameter; 
and  Sb  =  local  bottom  slope  given  by  Equation  11.  The  parameter  a  is  the 
ratio  between  the  wave  length  (2n/k)  and  the  horizontal  length  (3  h  /Sb) 
imposed  by  the  small  depth  and  relatively  steep  slope  where  the  lower  limit 
of  a  =  1  corresponds  to  the  formula  by  Battjes  and  Stive  (1985),  who  also 
assumed  Hb  =  Hm.  The  fraction  Q  is  zero  for  no  wave  breaking  and  unity 
when  all  waves  break.  The  requirement  of  o  <  Q  <  1  implies  Hrms  <  Hm  but 
Hrms  can  become  larger  than  H„,  in  very  shallow  water.  When  Hrms  >  Hm, 

Q  is  set  to  one  and  Hb  equal  to  Hrms.  In  addition,  the  upper  limit  of 
o,  =oj]/h  from  Equation  31  is  imposed  as  a,  <  1  in  very  shallow  water 

(Kobayashi  et  al.  1998).  The  breaker  ratio  parameter  y  in  Equation  38  is 
typically  in  the  range  of  y  =  0.5  -  1.0  (Kobayashi  et  al.  2007a)  but  should  be 
calibrated  to  obtain  a  good  agreement  with  the  measured  cross-shore 
variation  of  on  if  such  data  are  available.  The  default  value  for  field  applica¬ 
tions  is  y  =  0.5.  An  option  is  provided  in  CSHORE  to  estimate  y  with  the 
empirical  formula  developed  by  Apotsos  et  al.  (2008)  using  field  data. 

On  the  other  hand,  the  energy  dissipation  rate  Df  due  to  bottom  friction  in 
Equation  36  is  expressed  as: 


Df  =  |b fbU3.  (39) 

Substitution  of  Ua  given  in  Equation  27  into  Equation  39  yields: 

Df  =  \pfbo3TGf  ;Gf  =  J F3f(r)dr  (40) 

— oo 

where /(r)  is  given  by  Equation  29. 

The  energy  equation  for  the  roller  given  by  Equation  9  reduces  to  that  used 
by  Ruessink  et  al.  (2001)  for  the  case  of  alongshore  uniformity: 


^(pC2<?r  cos0)  =  Db  -  Dr  ;  Dr  =  pgfirqr 


(41) 
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where  the  roller  dissipation  rate  D,  is  assumed  to  equal  the  rate  of  work  to 
maintain  the  roller  on  the  wave-front  slope  (3r  which  is  of  the  order  of  o.i. 
Use  is  made  of  the  empirical  formula  given  by  Equation  10  proposed  by 
Kobayashi  et  al.  (2005)  who  included  the  local  bottom  slope  effect.  If  the 
roller  is  neglected,  qr=  o  and  Equation  41  yields  Dr  =  Db.  The  roller  effect 
improves  the  agreement  between  computed  and  observed  for  the  longshore 
current  (Kobayashi  et  al.  2007a). 

Equations  19-41  are  the  same  as  those  used  by  Kobayashi  et  al.  (2007a) 
who  assumed  Qx  =  q0  =  o  in  Equation  19  and  neglected  the  wind  shear 
stresses  in  Equations  22  and  23,  and  used  linear  shallow-water  wave 
theory  with  C  =  (gh  )°s  in  Equation  30.  Substitution  of  Equation  31  and 
32  into  Equation  19  yields: 


fr  9^ 

U  =  -^ouo, 


Cqr 


90, 


4 


(42) 


The  landward-marching  computation  starting  from  x  =  o  outside  the  surf 
zone  is  the  same  as  that  of  Kobayashi  et  al.  (2007a). 


Approximate  analytical  equations  of  Gbx,  Gby  and  G/ given  by  Equations  34 
and  40  are  obtained  by  Kobayashi  (2009a)  to  reduce  the  computation  time 
and  improve  the  numerical  stability.  The  function  Fa  given  in  Equation  27 
with  Equation  28  is  rewritten  as: 


F.  = 


+f: 


0.5 


(43) 


with 


rm  =  -  (U,  cos  9  +  K  sin  9) ;  Fm  =  V *  cos  9  -  U,  sin  9  (44) 

Equation  43  is  approximated  as: 

^a  =  (r-O  +  |4|forr>0; 

Fa=-{r-rm)  +  \Fm\^Vr<0 


(45) 
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Substituting  Equation  45  into  Equations  34  and  40  and  integrating  the 
resulting  equations  analytically,  we  obtain  approximate  expressions  for 
G bx,  Gby  and  Gf. 


GiI  =  ^|y.-r„,c0se)  +  t/,|F„1|;  (46) 

Gi„=]|(V,-r„sin0)  +  K.|F„,|;  (47) 

G,  =  2  J| +(l  +  t/,2  +V.2)|F„,|  +  +V.2  +  2r2)  (48) 

which  depends  on  sin0  (cos  0  >  o  assumed),  rm  and  Fm  where  Equation  44 
yields  U,  =  -  (rm  cos  0  +  Fm  sin  0)  and  V,  =  ( Fm  cos  0  -  rm  sin  0). 

For  the  case  of  normally  incident  waves  with  no  wind,  sin  0  =  o  and  U# 

=  o.  Equations  46-48  yield  Gbx  =  1.6  U*  ,  Gby  =  o,  and  G/  =  (1.6  +  2.4  JJ; ). 

Equation  23  requires  ny  =  o  for  Qx  =  o  (no  wave  overtopping)  and 
Equation  33  yields  Gby  =  o.  As  a  result,  Equation  47  is  exact.  For  sin  0  =  o 
and  V,  =  o,  Gbx  and  G/,  given  by  Equations  34  and  40,  can  be  integrated 

analytically  as  presented  by  Kobayashi  et  al.  (2007b)  who  approximated 
the  analytical  expressions  of  Gbx  and  G/as  Gbx  =  1.64  U,  and  G/=  (1.6  + 

2.6  U* ).  These  approximate  equations  are  very  similar  to  the  above 

equations  obtained  from  Equations  46  and  48.  For  the  case  of  |  sin  0|  <C  1 
and  |  G,  |  <SC  |  V,  |,  Equation  47  can  be  approximated  as  Gby  =  V,  (0.8  +  | 

V,  |).  Using  field  data  and  probabilistic  analyses,  Feddersen  et  al.  (2000) 
obtained  Gby  =  V,  (1.162  +  V2  )°-5.  The  difference  between  these  two 
approximate  equations  for  Gby  is  less  than  20  percent  for  |  V*  \  <  1.4,  which 
is  typically  satisfied. 

Kobayashi  et  al.  (2009a)  compared  the  approximate  values  of  Gbx,  Gby  and 
Gf  given  by  Equations  46-48  with  the  exact  values  of  Gbx,  Gby  and  G/ 
obtained  by  the  numerical  integration  of  Equations  34  and  40.  The 
percentage  error  was  typically  about  10  percent  and  always  less  than 
35  percent  for  the  ranges  of  |sin  0|  <  1,  |  rm  \  <  1  and  |Fm|  <  1.  This  error  is 
probably  less  than  the  uncertainty  of  the  bottom  friction  factor  fb. 
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5  Sediment  Transport  Model  in  Wet  Zone 


The  combined  wave  and  current  model  CSHORE  predicts  the  spatial 
variations  of  the  hydrodynamic  variables  used  in  the  following  sediment 
transport  model  for  given  beach  profile,  water  level  and  seaward  wave 
conditions  at  *  =  o.  The  bottom  sediment  is  assumed  to  be  uniform  and 
characterized  by  d50  =  median  diameter;  Wf  =  sediment  fall  velocity;  and 
s  =  sediment  specific  gravity.  The  sediment  transport  model  developed  for 
CSHORE  is  modified  slightly  for  the  horizontally  2-D  model  C2SHORE. 


The  spatial  variation  of  the  degree  of  sediment  movement  is  estimated 
using  the  critical  Shields  parameter  \pc  (Madsen  and  Grant  1976)  which  is 
taken  as  \j/c  =  0.05.  The  instantaneous  bottom  shear  stress  t  h  is  assumed 
to  be  given  by  x  i,  =  0.5  pfi,  U2a  with  U„  given  in  Equation  26.  The  sediment 
movement  is  assumed  to  occur  when  x  i,  exceeds  the  critical  shear  stress, 
pp(s  -i)d50  \|/c.  The  probability  Pb  of  sediment  movement  can  be  shown  to 
be  the  same  as  the  probability  of  (r  =  where 

Rb  =  [2  g  (s-t)  d5 o  \|/c/6"1]0-5/ar  and  rm  and  Fm  are  defined  in  Equation  44. 
For  the  Gaussian  variable  r  given  by  Equation  29,  Pb  is  given  by: 


pb  =  ~erfc 


Fb~r,r 


+—erfc 

2 


Fh+r„ 


for  F?  >  0 


(49) 


and  Pb  =  1  for  <  o  where  erfc  is  the  complementary  error  function.  The 

value  of  Pb  computed  from  x  =  o  located  outside  the  surf  zone  increases 
landward  and  fluctuates  in  the  surf  and  swash  zones,  depending  on  the 
presence  of  a  bar  or  a  terrace  that  increases  the  local  fluid  velocity. 

The  spatial  variation  of  the  degree  of  sediment  suspension  is  estimated 
using  the  experimental  finding  of  Kobayashi  et  al.  (2005)  who  showed  that 
the  turbulent  velocities  measured  in  the  vicinity  of  the  bottom  were  related 
to  the  energy  dissipation  rate  due  to  bottom  friction.  Representing  the 
magnitude  of  the  instantaneous  turbulent  velocity  by  (D'f/ p)J/3  with 
D’f  =  0.5  pfb  U3  in  light  of  Equation  39,  the  probability  Ps  of  sediment 

suspension  is  assumed  to  be  the  same  as  the  probability  of  (D'f/p)1/ 3 
exceeding  the  sediment  fall  velocity  w/.  The  probability  Ps  is  then  equal  to 
the  probability  of  F3  =  {R't,  -  F3 )  with  Rs  =  [(2//),) ‘/s  Wf/ori  and  is  given  by: 
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Ps  =  —erfc 

S  2 


F  —r 

s  m 


1  r 

+— ertc 

2 


F  +r 

s  1  m 


for  Ff  >  0 


(50) 


and  Ps  =  l  for  Fs2  <  o.  If  Ps  >  Pb,  use  is  made  of  Ps  =  Pb  assuming  that 

sediment  suspension  occurs  only  when  sediment  movement  occurs.  Fine 
sands  on  beaches  tend  to  be  suspended  once  their  movement  is  initiated. 


The  suspended  sediment  volume  Vs  per  unit  horizontal  bottom  area  is 
estimated  by  modifying  the  sediment  suspension  model  by  Kobayashi  and 
Johnson  (2001): 


v  _p  eBDr+efDf 

Spg(s-l)wf 


„  \0.5  /  n  \  0.5  dz,  dz 

l  +  SL  1  +  S£  ;Sta=^;s,,„  = 


dx 


dy 


(51) 


where  S bx  =  cross-shore  bottom  slope;  Shy  =  longshore  bottom  slope;  and  ep 
and  e/  =  suspension  efficiencies  for  the  energy  dissipation  rates  D,  and  Df 
due  to  wave  breaking  and  bottom  friction,  respectively.  Use  is  made  of 
es  =  0.005  and  ef~  0.01  as  typical  values  in  the  computation  of  berm  and 
dune  erosion  but  the  value  of  ep,  is  uncertain  and  should  be  calibrated  if  Vs  is 
measured  (Kobayashi  et  al.  2007a).  The  sediment  suspension  probability  Ps 
is  added  to  Equation  51  to  ensure  Vs  =  o  if  Ps  =  o.  The  term  involving  Sbx  and 
Sby  is  the  actual  bottom  area  per  unit  horizontal  bottom  area  and  essentially 
unity  except  for  very  steep  slopes.  For  the  case  of  alongshore  uniformity, 

Sby  =  o.  The  cross-shore  and  longshore  suspended  sediment  transport  rates 
qsx  and  qsy  are  expressed  as: 


Qsx  =  axUVs ;  -  W„ ;  a,.  = 


lsy 


a  +  O^/tancp)0-5 


>  a 


(52) 


where  a  =  empirical  suspended  load  parameter  and  cj)  =  angle  of  internal 
friction  of  the  sediment  with  tan(j)  =  0.63  for  sand  (Bailard  1981).  The 
parameter  a  accounts  for  the  onshore  suspended  sediment  transport  due  to 
the  positive  correlation  between  the  time-varying  cross-shore  velocity  and 
suspended  sediment  concentration.  The  value  of  a  increases  to  unity  as  the 
positive  correlation  decreases  to  zero.  For  the  three  small-scale  equilibrium 
profile  tests  conducted  by  Kobayashi  et  al.  (2005),  a  was  of  the  order  of  0.2. 
The  effect  of  the  cross-shore  bottom  slope  on  ax  was  included  by  Kobayashi 
et  al.  (2009b)  to  increase  berm  and  dune  erosion.  For  Sbx  <  o ,ax  =  a.  The 
cross-shore  suspended  sediment  transport  rate  qsx  is  negative  (offshore) 
because  the  return  (undertow)  current  U  is  negative  (offshore).  On  the 
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other  hand,  the  longshore  suspended  sediment  transport  rate  qsy  in 
Equation  52  neglects  the  correlation  between  the  time-varying  longshore 
velocity  and  suspended  sediment  concentration,  which  appears  to  be  very 
small  if  the  longshore  current  V  is  sufficiently  large.  Payo  et  al.  (2009) 
verified  Equation  52  using  velocities  and  sand  concentrations  measured 
along  20  transects  at  the  Field  Research  Facility  (FRF)  at  Duck,  North 
Carolina  during  a  storm  in  1997. 

The  formulas  for  the  cross-shore  and  longshore  bedload  transport  rates 
qbx  and  qby  are  devised  somewhat  intuitively  because  bedload  in  the  surf 

zone  has  never  been  measured.  The  time-averaged  rates  qbx  and  qblJ  are 

tentatively  expressed  as: 

qtx  =  BlJF+V^ju ;  qby  =  J?(tF Tv^jv  (53) 


where  B  =  empirical  parameter.  Equation  53  may  be  regarded  as  a  quasi¬ 
steady  application  of  the  formula  of  Meyer-Peter  and  Mueller  (e.g., 
Ribberink  1998).  Substitution  of  U  and  V given  in  Equation  27  with 
Equations  28  and  29  into  Equation  53  yields: 


qbx  =  Bo3T[h  +  £/,K2  +  2  Fm  sinG 


9by  =  Bai 


VJl  +  U?  +  K2 ) - 2r  sinG 


(54) 

(55) 


where  b,  =  (3U,  + 1/,3)  and  Fm  and  rm  are  defined  in  Equation  44. 


Equation  54  and  55  yield  qbx  =  KBo3T  and  qhy  =  o  for  normally  incident 
waves  with  sin  9  =  o  and  K  =  o.  The  expressions  of  B  and  are  obtained 

by  requiring  that  qbx  =  KBo3T  reduces  to  the  onshore  bedload  formula 

proposed  by  Kobayashi  et  al.  (2008a)  for  normally  incident  waves,  which 
synthesized  existing  data.  The  proposed  formulas  are  written  as: 

bP 

Qto  =-r^7-yTb  +  U.V.::+2F„sme)G,(Sb)  (56) 

9{s  -Lj 


9bu  ~ 


bPh 


9[s- 1 


a: 


V, (l  +  U?+V.2  -2r  sinG 


GAS ; 


Jby 


(57) 
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where  b  =  empirical  bedload  parameter;  and  Gs  =  bottom  slope  function. 
The  sediment  movement  probability  Pb  given  in  Equation  49  accounts  for 
the  initiation  of  sediment  movement.  It  is  noted  that  6*  =  1  in  Equation  56 

to  compensate  for  the  limitations  of  Equation  53  and  the  Gaussian  distribu¬ 
tion  of  the  horizontal  velocity  used  in  Equations  28  and  29,  as  discussed  by 
Kobayashi  et  al.  (2008a).  They  calibrated  b  =  0.002  using  the  20  water 
tunnel  tests  of  Ribberink  and  Al-Salem  (1994),  the  four  large-scale  wave 
flume  tests  of  Dohmen-Janssen  and  Hanes  (2002),  and  the  24  sheet  flow 
tests  by  Dohmen-Janssen  et  al.  (2002).  Furthermore,  this  simple  bedload 
formula  is  consistent  with  the  sheet  flow  model  for  onshore  bar  migration 
by  Trowbridge  and  Young  (1989)  and  the  energetics-based  bedload  formula 
for  steady  flow  by  Bagnold  (1966)  if  the  steady  flow  formula  is  applied  in  the 
time-averaged  manner.  The  onshore  bedload  transport  predicted  by 
Equation  56  is  consistent  with  the  field  observations  of  onshore  ripple 
migration  by  Becker  et  al.  (2007)  and  Masselink  et  al.  (2007).  The  offshore 
suspended  sediment  transport  predicted  by  Equation  52  is  consistent  with 
the  field  measurement  during  a  storm  by  Madsen  et  al.  (1994).  The  condi¬ 
tion  of  ( qbx  +  qsx )  =  0  for  an  equilibrium  profile  along  with  additional 

assumptions  can  be  shown  to  yield  the  equilibrium  profile  popularized  by 
Dean  (1991). 

The  bottom  slope  function  Gs(Shx)  was  introduced  by  Kobayashi  et  al. 
(2008a)  to  account  for  the  effect  of  the  steep  cross-shore  slope  Sbx  on  the 
bedload  transport  rate  and  is  expressed  as: 

Gs  ( Sbx )  =  tan  (p  /  (tan  (p  +  Sbx )  for  -  tancp  <  Sbx  <  0  (58) 

Gs  (sta )  =  (tan cp  - 2Ste )  /  (tan  cp  -  ) for  o<Sbx<  tancp  (59) 

where  Gs  =  1  for  Sbx  =  o.  Equation  58  corresponds  to  the  functional  form  of 
Gs  used  by  Bagnold  (1966)  for  steady  stream  flow  on  a  downward  slope 
with  Sbx  <  o  where  the  downward  slope  increases  qbx.  Equation  59  ensures 
that  Gs  approaches  negative  infinity  as  the  upward  slope  Sbx  approaches 
tan  (j).  Equations  58  and  59  reduce  to  Gs  =  (1  -  SW/tan  (j))  for  \Sbx\  <  tan  cj). 
Equation  56  with  Gs  given  by  Equations  58  and  59  implies  that  the 
bedload  transport  rate  qbx  is  positive  (onshore)  for  Sbx  <  (tan  §)/ 2.  and 

negative  (offshore)  for  Sbx  >  (tan  c)))/2.  Use  is  made  of  |  Gs\  <  Gm  =  10  to 
avoid  an  infinite  value  in  the  computation.  The  computed  profile  change  is 
not  very  sensitive  to  the  assumed  value  of  Gm  because  the  beach  profile 
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changes  in  such  a  way  to  reduce  a  very  steep  slope  except  in  the  region  of 
scarping  (e.g.,  Seymour  et  al.  2005).  The  effect  of  the  longshore  bottom 
slope  Sby  is  included  in  Equation  57  using  the  same  bottom  slope  function 
Gs(Sby )  but  has  never  been  validated  for  lack  of  suitable  data. 

The  landward  marching  computation  of  the  present  time-averaged  model 
ends  at  the  cross-shore  location  x  =  xm  where  the  mean  water  depth  h  is 
less  than  1  cm.  No  reliable  data  exists  for  suspended  sand  and  bedload 
transport  rates  in  the  zone  which  is  wet  and  dry  intermittently.  Conse¬ 
quently,  the  following  simple  procedure  was  proposed  by  Kobayashi  et  al. 
(2008a)  to  deal  with  the  zone  with  a  bottom  slope  Sbx  >  tan  4>.  The  cross¬ 
shore  total  sediment  transport  rate  qx  =  (qSx  +  qbx )  at  x  =  xm  is  denoted  by 
Qxm»  If  q  xm  is  negative  (offshore),  qx  is  extrapolated  linearly  to  estimate  qx  on 
the  scarped  face  with  Sbx  >  tan  (j)  such  that: 

Qx  =  Qxxn  (*e  - x) / {xe  —  xm )  for  xm<x<  xe  (60) 

where  xe  =  landward  limit  of  the  scarping  zone  with  Sbx  >  tan  (j).  The 
extrapolated  qx  is  in  the  range  of  qxm  <  qx  <  o  and  the  scarping  zone  is 

eroded  due  to  the  offshore  sediment  transport.  This  simple  procedure  is 
effective  for  a  high  and  wide  dune,  that  is  typical  in  the  Netherlands  (e.g., 
van  Gent  et  al.  2006),  but  does  not  allow  onshore  sediment  transport  due 
to  overwash.  As  a  result,  no  wave  overtopping  has  been  allowed  so  far  and 
q0  =  o  in  Equation  19. 

The  beach  profile  change  is  computed  using  the  continuity  equation  of 
bottom  sediment: 


l~np)— 

p)  dt 


Ox  dy 


(61) 


where  np  =  porosity  of  the  bottom  sediment  which  is  normally  taken  as 
np  =  0.4;  t  =  slow  morphological  time  for  the  change  of  the  bottom 
elevation  Zb;  and  qy  =  ( qsy  +  qby)  =  longshore  total  sediment  transport  rate. 
For  the  case  of  alongshore  uniformity,  the  third  term  in  Equation  61  is 
zero.  Equation  61  is  solved  using  an  explicit  Lax-Wendroff  numerical 
scheme  (e.g.,  Nairn  and  Southgate  1993)  to  obtain  the  bottom  elevation  at 
the  next  time  level.  This  computation  procedure  is  repeated  starting  from 
the  initial  bottom  profile  until  the  end  of  a  profile  evolution  test. 
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6  Permeable  Layer  Model  in  Wet  Zone 

The  combined  wave  and  current  model  CSHORE  is  extended  to  allow  the 
presence  of  a  permeable  layer  in  the  computation  domain.  Figure  2  shows 
an  example  of  irregular  wave  overtopping  of  a  permeable  slope  where  x  = 
onshore  coordinate:  z  =  vertical  coordinate,  r|  =  mean  free  surface  elevation 
above  SWL;  Mean  Water  Level  (MWL)  is  the  SWL  +  q ;  S  =  storm  tide 

above  z  =  o;zb  =  bottom  elevation;  h  =  mean  water  depth;  U  =  instant- 
neous  depth-averaged  cross-shore  velocity  above  the  bottom;  zp  =  elevation 
of  the  lower  boundary  of  the  permeable  layer;  hp  =  (zb  -zp)  =  vertical 
thickness  of  the  permeable  layer;  and  Up  =  instantaneous  cross-shore 
discharge  velocity  inside  the  permeable  layer.  The  cross-shore  profiles  of 
z/,(x)  and  zp(x)  are  specified  as  input  where  hp  =  o  in  the  zone  of  no 
permeable  layer.  The  lower  boundary  located  at  z  =  zp  is  assumed  to  be 
impermeable  for  simplicity.  Kobayashi  et  al.  (2007b)  developed  a 
permeable  layer  model  in  the  wet  zone  for  normally  incident  waves.  This 
model  is  extended  to  obliquely  incident  waves  in  the  following  but  the 
extended  model  has  not  yet  been  verified. 


The  time-dependent  model  for  the  flow  over  a  permeable  layer  in  shallow 
water  developed  by  Kobayashi  and  Wurjanto  (1990)  and  Wurjanto  and 
Kobayashi  (1993)  is  time-averaged  and  simplified  to  account  for  the 
permeable  layer  in  the  cross-shore  model  CSHORE.  The  vertically- 
integrated  continuity  Equation  19  is  modified  as: 
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Qx  =  hU  +  ^—^cosB  +  qrcosQ;Qx  +  hpUp  =  q0  (62) 


where  Up  =  time-averaged  cross-shore  discharge  velocity;  [hpUpj  =  water 

flux  inside  the  permeable  layer  with  its  vertical  thickness  hp  ;  and  qa  = 

combined  wave  overtopping  rate  above  and  through  the  permeable  layer. 
The  cross-shore  and  longshore  momentum  Equations  22  and  23  are 
assumed  to  remain  the  same,  neglecting  the  momentum  fluxes  into  and 
out  of  the  permeable  layer  where  the  bottom  friction  factor  fb  for  xbx  and 

xhy  given  by  Equation  33  includes  the  effect  of  the  surface  roughness  of 
the  permeable  layer  and  was  calibrated  in  the  range  of  fb  =  0.01  -  0.05 

(Kobayashi  et  al.  2007b).  For  the  case  of  alongshore  uniformity  and 
negligible  momentum  fluxes  into  and  out  of  the  permeable  layer,  the  time- 
averaged  longshore  discharge  velocity  Vp  is  assumed  to  be  zero  because  of 

no,  or  negligible,  driving  force  to  cause  the  longshore  discharge  inside  the 
permeable  layer.  It  is  noted  that  the  assumption  of  V  =  0  cannot  be 

validated  at  present  for  lack  of  suitable  data. 


On  the  other  hand,  the  wave  action  Equation  36  is  modified  as: 


d_ 

dx 


E_ 

to 


Cffcos9  + 


D B  +  Dj  +  Dp 

co 


(63) 


where  D  =  energy  dissipation  rate  due  to  flow  resistance  in  the  permeable 

layer,  assuming  that  the  energy  influx  into  the  permeable  layer  equals  the 
dissipation  rate  Dp  per  unit  horizontal  area.  The  dissipation  rate  Dp  is 

expressed  as  (Wurjanto  and  Kobayashi  1993): 


Dv  =  PK 


a  I Uz+Vz 

p\  p  p 


+  P 


<UZ+VZ 

P\  P  '  P 


1.5 


(64) 


where  a  and  pp  =  laminar  and  turbulent  flow  resistance  coefficients, 
respectively,  and  V  =  instantaneous  longshore  discharge  velocity. 
Kobayashi  et  al.  (2007b)  modified  the  formulas  for  a  and  pp  proposed 
by  van  Gent  (1995)  as  follows: 
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with 


S  =  ao 


n. 


v 


P2 


n. 


Dz  A=e.+f 

^nSO  Wp 


Pi 


?oH. R  _7-5Po(1-np) 
n^n5o  ’  2  V2njT 


(65) 


(66) 


where  a0  and  P0  =  empirical  parameters  calibrated  as  a0  =  1,000  and  P0 
=  5 ;  np  =  porosity  of  the  permeable  layer  consisting  of  stone;  v  = 
kinematic  viscosity  of  the  fluid;  Dn 50  =  nominal  stone  diameter  defined  as 
Dn50  =  {Mso/Psf/3  with  M5o  =  median  stone  mass  and  ps  =  stone  density; 
op  =  standard  deviation  of  the  instantaneous  discharge  velocity;  and 
T  =  intrinsic  wave  period  used  in  Equation  38. 

The  discharge  velocities  Up  and  Vp  in  Equation  64  are  assumed  to  be 
expressed  as: 


Up  =  Up  +ro p  cosG ;  Vp  =  rap  sinG;  (67) 

where  r  =  Gaussian  variable  whose  probability  density  function  is  given  by 
Equation  29;  and  0  =  incident  wave  angle  for  the  oscillatory  velocity 
direction  above  and  inside  the  permeable  layer.  The  assumptions  of  the 
Gaussian  velocity  distribution  and  V  =  0  allow  one  to  represent  the 

discharge  velocities  by  the  mean  cross-shore  discharge  velocity  U  and  the 
standard  deviation  op .  Substitution  of  Equation  67  into  Equation  64  yields: 


1 - 2  0  \  [2  /  \ 

2  1 

dp=pk- 

ap[uP  +ap)  +  ^-(p2  +  P1ap) 

2o2p+Up  (l  +  2cos20) 

where  use  is  made  of  the  approximate  expression  of  Gf  given  by 


Equation  48  and  the  assumption  of 


£/  sin0 


<g?  rs 


tn  cirrmlifu  'Rmintinn 


Approximate  equations  for  Up  and  a  are  derived  in  the  following. 
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Neglecting  the  inertia  terms  in  the  cross-shore  momentum  equation  for 
the  flow  inside  the  permeable  layer  (Kobayashi  and  Wurjanto  1990),  the 
local  force  balance  between  the  cross-shore  hydrostatic  pressure  gradient 
and  flow  resistance  is  assumed  to  be: 

gfx  +  apUM{Ul+Vir=0  (69) 

Equation  69  is  averaged  probabilistically  using  Equation  67.  For  the  case  of 
alongshore  uniformity,  the  averaged  force  balance  equation  is  expressed  as: 

9fa+Up  aP+)/|'(P2+Piap)(1  +  COs20)  =0  (7Q) 

where  use  is  made  of  the  approximate  expression  of  Gbx  given  by 
Equation  46  and  the  assumption  of  Up  sin0  <C  op  to  simplify  Equation  70. 

It  is  noted  that  the  local  force  balance  between  the  longshore  hydrostatic 
pressure  gradient  and  flow  resistance  yields  V  =  0  for  the  case  of  along¬ 
shore  uniformity  where  q  is  independent  of  the  longshore  coordinate  y .  To 
derive  an  equation  for  a  ,  the  approximate  analytical  method  used  by 
Kobayashi  et  al.  (2007b)  is  adopted.  Equation  69  is  linearized  as: 

+K+1-9pp°pK=0  ™ 

which  is  used  to  obtain: 

ap+1.9(p2+$1op)}op=gkho.;o.=o1i/h  (72) 

where  the  wave  number  k  is  computed  using  Equation  2.  Equation  72  can 
be  solved  analytically  to  obtain  op  for  known  kho , .  After  op  is  obtained, 

Equation  70  is  used  to  calculate  Up  for  known  dr\/dx.  The  energy  dissipa¬ 
tion  rate  Dp  is  computed  using  Equation  68.  Equation  62  for  assumed  qo  is 

used  to  obtain  Qx  and  U  where  U  is  expressed  by  Equation  42. 
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7  Irregular  Wave  Runup  and  Overtopping 

The  time-averaged  model  CSHORE  does  not  predict  the  shoreline  oscillations 
on  beaches  and  coastal  structures  unlike  time-dependent  models  (e.g., 
Wuijanto  and  Kobayashi  1993).  To  compensate  for  this  shortcoming  of 
CSHORE,  Kobayashi  et  al.  (2008b)  proposed  a  probabilistic  model  for 
irregular  wave  runup  as  illustrated  in  Figure  3.  The  shoreline  oscillation  is 
assumed  to  be  measured  by  a  runup  wire  (RW  in  Figure  3)  placed  parallel  to 
the  bottom  elevation  z/,  at  a  vertical  height  of  Sr.  The  runup  wire  measures  the 
instantaneous  elevation  //,■  above  SWL  of  the  intersection  between  the  wire 
and  the  free  surface  elevation.  The  mean  rp  and  standard  deviation  ar  of  nr 

are  estimated  using  the  computed  cross-shore  variations  of  r|(x)  and  op  (x) 
of  the  free  surface  elevation  r\  above  SWL.  The  probabilities  of  //,-  exceeding 
(r\r  +  or  j ,  rp ,  and  (rp  -  a,,  j  are  assumed  to  be  the  same  as  the  probabilities 

of  r|  exceeding  (13  +  aJ ,  13 ,  and  (13  —  a n  j ,  respectively.  The  elevations  of  Zh 

Z2,  and  Z3  of  the  intersections  of  (r\  +  an  ) ,  13 ,  and  (13  —  an )  with  the  runup 

wire  are  obtained  for  the  given  wire  elevation  (z/,  +  Sr).  The  obtained 
elevations  are  assumed  to  correspond  to  '  >lr+0r  !  =  Hr  ,  and 

Z3  =  ( rp  —  or  j .  The  mean  and  standard  deviation  of  13,.  are  estimated  as: 

fp  =  (Z1  +  Z2  +  Z3 )  /  3  ;  or  -  (Z1  -  Z3 )  /  2  (73) 

In  CSHORE,  13  and  on  are  replaced  by  (P10/2  +  zfaj  and  Pwo n  for  the 

computation  of  Zh  Z2  and  Z3  to  account  for  the  transition  from  the  wet 
zone  ( Pw  =  1 )  to  the  wet  and  dry  zone  ( Pw  <  1)  where  Pw  is  the  wet 
probability  explained  in  Chapter  8. 

The  runup  height  R  is  defined  as  the  crest  height  above  SWL  of  the 
temporal  variation  of  rp-  The  probability  distribution  of  linear  wave  crests 
above  the  mean  water  level  (MWL)  is  normally  given  by  the  Rayleigh 
distribution.  For  the  case  of  no  wave  overtopping,  the  runup  height 
IK  —  rp  j  above  the  mean  elevation  rp  is  assumed  to  be  given  by  the 

Rayleigh  distribution  (Kobayashi  et  al.  2008b).  The  exceedance 
probability  P(R)  of  the  runup  height  R  above  SWL  is  given  as: 
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Figure  3.  Definition  sketch  for  probabilistic  model  for  irregular  wave  runup. 


P(R)  =  exp 


-2 


R-^r 


ri,3  -a 


(74) 


where  i?i/3  =  significant  runup  height  which  is  defined  as  the  average  of  1/3 
highest  values  of  R. 

The  mean  r\r  related  to  wave  setup  is  normally  neglected  in  Equation  74 

for  the  prediction  of  irregular  wave  runup  on  steep  coastal  structures.  For 
the  1/5  and  1/2  permeable  slope  experiments  conducted  by  Kobayashi 
et  al.  (2008b),  Ri/3  was  estimated  as: 

^/3  =hr+(2  +  tan0)ar  (75) 

where  9  =  seaward  slope  angle  from  the  horizontal  and  tan  0  =  0.2  and  0.5 
in  the  experiments.  It  is  cautioned  that  Equations  74  and  75  have  been 
calibrated  only  for  permeable  slopes  with  tan  9  =  0.2  -  0.5  in  the  absence 
of  wave  overtopping. 

Wave  overtopping  occurs  when  the  individual  runup  height  R  above  SWL 
exceeds  the  structure  crest  height  Rc  above  SWL  as  depicted  in  Figure  3. 
Wave  overtopping  reduces  R  exceeding  Rc  because  of  overtopping  flow  on 
the  crest.  Kobayashi  and  de  los  Santos  (2007)  adopted  the  following 
Weibull  distribution: 
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/  —  \K 

P(R)  =  ex  p-2  ^~n!L  (76) 

["L/3  Hr  J 

with 

K  =  2  +  0.5ft-3  ;  R,  =  (R,  -  n, ) /  (Rll3  -  nr)  (77) 

where  k  =  shape  parameter  with  k  =  2  for  the  Rayleigh  distribution  given 
by  Equation  74;  and  R  =  normalized  crest  height  related  to  the  wave 

overtopping  probability  P0.  The  probability  P0  of  R  exceeding  R(  in 
Equation  76  is  given  by: 

P0=exp(-2P*K)  (78) 

It  should  be  noted  that  the  empirical  formula  for  k  given  by  Equation  77 
has  been  calibrated  using  only  22  permeable  slope  tests  so  far.  The 
formula  for  R,/3  given  by  Equation  75  has  been  found  to  be  applicable  to 
these  22  tests.  The  runup  height  R2%  for  the  two  percent  exceedance 
probability  obtained  using  Equation  76  is  given  by 

7?»  =  nr+(1.40f77?1/3-nr)  (79) 

where  the  shape  parameter  k  given  by  Equation  77  accounts  for  the 
decrease  of  R2%  due  to  the  decrease  of  the  normalized  crest  height  P,  and 

the  resulting  increase  of  the  wave  overtopping  probability  P0  given  by 
Equation  78. 

The  wave  overtopping  rate  q0  in  Equation  19  for  an  impermeable  slope  and 
in  Equation  62  for  a  permeable  slope  needs  to  be  estimated  if  wave 
overtopping  occurs  at  the  landward  end  of  the  computation  domain 
located  at  *  =  xe  in  Figure  2.  For  permeable  slopes,  Kobayashi  and 
de  los  Santos  (2007)  proposed  the  following  empirical  formula: 

Q0=^(P0)b'  qSWL+qs  (80) 

with 

go2 

Qswl  =  -^cosG  at  X  =  XSWL  (81) 
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where  a *  and  b,  =  empirical  parameters;  P0  =  wave  overtopping 
probability;  qSWL  =  wave-induced  onshore  flux  in  Equation  62  evaluated  at 
the  still  water  shoreline  located  at  x  =  xSWL  with  zb  (xSWL)  =  S  in  Figure  2; 
and  qs  =  seepage  rate  through  the  permeable  layer.  It  is  noted  that  the 

roller  effect  has  been  neglected  for  permeable  slopes  because  of  its 
negligible  effect  and  qr  =  0  in  Equation  62.  The  empirical  parameters  a* 

and  b  are  assumed  to  depend  on  the  horizontal  width  Lh  of  the  permeable 
surface  above  the  upper  limit  of  wave  setup  located  at  (xr ,  zr )  in  Figure  3, 

where  the  infiltration  of  overtopped  water  is  assumed  to  be  vertical  due  to 
gravity.  The  empirical  formulas  based  on  32  tests  were  expressed  as: 

a*exp(-0.1L*);&  =  1  +  0.11.  ;L»  =  Lh  / Dn50  (82) 

where  L*  =  infiltration  width  normalized  by  the  nominal  stone  diameter 

Dn5o,  crudely  representing  the  horizontal  number  of  stones  above  the 
maximum  wave  setup. 

On  the  other  hand,  Kobayashi  and  de  los  Santos  (2007)  estimated  the 
seepage  rate  qs  for  normally  incident  waves  as: 


Qs 


0.2(zr 


for  zr  >  zf, 


(83) 


where  ze  =  elevation  of  the  landward  end  of  the  impermeable  surface  z 
as  shown  in  Figure  2;  and  p1  =  turbulent  flow  resistance  coefficient 

defined  in  Equation  66.  To  derive  Equation  83,  the  seepage  flow  was 
assumed  to  be  driven  by  the  horizontal  pressure  gradient  from  the  point 
(xr,zr)  to  the  point  ( xe ,  ze ) .  Consequently,  qs  =  0  if  zr  <  ze .  If  xr=xe,  the 

permeable  layer  is  always  wet  and  qs  =hpUP  at  x  =  xe  where  the  water 
flux  hpUP  in  the  permeable  layer  is  included  in  the  continuity  equation 
given  in  Equation  62. 

Kobayashi  et  al.  (2007c)  examined  the  transition  from  little  wave 
overtopping  to  excessive  wave  overtopping  and  overflow  on  an  impermeable 
smooth  levee  with  a  seaward  slope  of  1/5  in  wave-flume  experiments 
consisting  of  107  tests.  For  the  impermeable  slope,  Equations  75  and  77 
used  for  the  permeable  slope  had  to  be  modified  as: 
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R1/3  =  nr+4ar  ;k  =  2  (84) 

The  wave  overtopping  probability  Po  is  given  by  Equation  78  with  k  =  2 
where  the  normalized  crest  height  R  above  SWL  is  defined  in 
Equation  77  with  R1/3  given  by  Equation  84.  It  is  noted  P  =  1  if  R  <  0  . 
For  the  impermeable  slope,  the  seepage  rate  qs  =  0  in  Equation  80  and 
Equation  82  yields  a*  =  1  and  b,  =  1  for  Lh  =  0 .  As  a  result,  the  wave 
overtopping  rate  qo  is  given  by  q0  =  PoqSWL  •  F°r  the  case  of  combined 

wave  overtopping  and  overflow,  Kobayashi  et  al.  (2007c)  expressed  the 
combined  rate  qa  as: 

Qo  =  PoQswl  d"  Hswl  gHSWL  for  HSWL  >  0  (85) 


with 


Hswl  ~  h  at  x  ~  xswl  (86) 

where  HSWL  =  head  for  the  overflow;  q  =  MWL  above  SWL;  and  R(  = 
levee  crest  height  above  SWL.  If  Rc  <  0 ,  the  levee  crest  is  below  SWL  and 
xSWL  is  chosen  at  the  seaward  edge  of  the  levee  crest.  For  HSWL  >  0 ,  HSWI 
is  the  MWL  above  the  levee  crest  and  gH SWL  may  be  regarded  as  the 
water  velocity  on  the  crest. 

In  summary,  Equations  73-86  are  essentially  empirical  and  are  used  in  the 
cross-shore  model  CSHORE  to  predict  irregular  wave  runup,  overtopping, 
seepage,  and  overflow  on  permeable  and  impermeable  structures.  These 
equations  have  not  been  verified  for  irregular  wave  overtopping  and 
overflow  of  dunes.  These  equations  do  not  predict  the  spatial  variations  of 
the  hydrodynamic  variables  required  for  the  sediment  model  and  the 
computation  of  dune  profile  evolution.  Consequently,  a  hydrodynamic 
model  for  the  intermittently  wet  zone  landward  of  the  maximum  wave  setup 
has  been  developed  in  Chapter  8.  The  formula  for  P0  given  by  Equation  78 
and  the  formulas  for  q0  given  by  Equations  80  and  85  are  removed  and 
replaced  by  new  formulas  for  P0  and  q0  based  on  the  hydrodynamic  model 
for  the  wet  and  dry  zone.  The  values  of  q ,  or ,  Rq3  and  R2%  are  computed  in 

CSHORE  because  this  hydrodynamic  model  does  not  predict  individual 
wave  runup  events. 
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8  Model  for  Impermeable  Wet  and  Dry  Zone 


Time-dependent  numerical  models  such  as  the  nonlinear  shallow-water 
wave  model  by  Kobayashi  et  al.  (1989)  can  predict  the  water  depth  and 
horizontal  velocity  in  the  intermittently  wet  and  dry  (swash)  zone  on 
beaches  and  inclined  structures.  However,  the  time-dependent  hydro- 
dynamic  computation  requires  considerable  computation  time  and  may  not 
lead  to  an  accurate  prediction  of  dune  profile  evolution  in  view  of  the  earlier 
attempt  by  Tega  and  Kobayashi  (1996).  A  time-averaged  probabilistic  model 
is  developed  here  to  predict  the  cross-shore  variations  of  the  wet  probability 
and  the  mean  and  standard  deviation  of  the  water  depth  and  cross-shore 
velocity  in  the  swash.  The  developed  model  is  very  efficient  computationally 
and  calibrated  using  a  large  number  of  data  sets.  The  present  model  is 
limited  to  normally  incident  waves  and  alongshore  uniformity.  A  sediment 
transport  model  in  the  swash  zone  is  formulated  by  modifying  the  sediment 
transport  model  in  the  wet  zone. 

Water  depth  and  velocity 

Van  Gent  (2002a)  and  Schiittrumpf  and  Oumeraci  (2005)  analyzed  the 
water  depth  and  velocity  of  waves  overtopping  of  dikes.  Kobayashi  et  al. 
(2010a)  expanded  their  analyses  for  the  prediction  of  wave  overtopping 
and  overwash  as  presented  in  the  following. 

For  normally  incident  waves  on  impermeable  beaches  and  inclined 
structures  of  alongshore  uniformity,  the  time-averaged  cross-shore 
continuity  and  momentum  equations  derived  from  the  nonlinear  shallow- 
water  wave  equations  are  expressed  as: 

hU  =  qo  (87) 


d_ 

dx 


hU2  +—h2 


-9Sbxh 


|/tpp;Sb 


dzb 

dx 


(88) 


where  h  and  U  =  instantaneous  water  depth  and  cross-shore  velocity, 
respectively;  qo  =  combined  wave  overtopping  and  overflow  rate;  g  =  gravi¬ 
tational  acceleration;  Sbx  =  cross-shore  bottom  slope;  and  fb  =  bottom 
friction  factor,  which  is  allowed  to  vary  spatially.  The  wave  energy  equation 
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corresponding  to  Equations  87  and  88  was  given  by  Kobayashi  and 
Wurjanto  (1992),  who  used  it  to  estimate  the  rate  of  wave  energy  dissipation 
due  to  wave  breaking.  The  wave  energy  equation  is  not  used  in  CSHORE 
because  no  formula  is  available  to  estimate  the  time-averaged  energy 
dissipation  rate  in  the  wet  and  dry  zone. 

The  instantaneous  water  depth  h  depends  on  the  cross-shore  coordinate  x. 
The  water  depth  h  at  given  x  is  described  probabilistically  rather  than  in  the 
time  domain.  Kobayashi  et  al.  (1998)  analyzed  the  probability  distributions 
of  the  free  surface  elevations  measured  in  the  shoaling,  surf  and  swash 
zones.  The  measured  probability  distributions  were  shown  to  be  in  agree¬ 
ment  with  the  exponential  gamma  distribution  which  reduces  to  the 
Gaussian  distribution  and  the  exponential  distribution  when  the  skewness 
approaches  zero  offshore  and  two  in  the  swash  zone,  respectively.  The 
assumption  for  the  Gaussian  distribution  assumed  in  Equation  29  has 
simplified  the  cross-shore  model  CSHORE  in  the  wet  zone  significantly.  The 
assumption  of  the  exponential  distribution  is  made  here  to  simplify  the 
cross-shore  model  in  the  wet  and  dry  zone.  The  probability  density  function 
f[h)  is  expressed  as: 


/(/i)iexp 


~p  i 

wh 


for  h  >  0 


with 


OO  OO 

Pw=  ff(h)dh;h  =  Jhf(h)dh 

0  0 


(89) 


(90) 


where  Pw  =  wet  probability  for  the  water  depth  h  >  o;  and  h  =  mean  water 
depth  for  the  wet  duration.  The  dry  probability  of  h  -  o  is  equal  to  (l  -  Pw ) 

.  The  mean  water  depth  for  the  entire  duration  is  equal  to  Pwh  .  The 
overbar  in  Equations  87  and  88  indicates  averaging  for  the  wet  duration 
only.  The  free  surface  elevation  (q  —  qj  above  MWL  is  equal  to  (h  —  h).  The 

standard  deviations  of  q  and  h  are  the  same  and  given  by 


o 


n 


h 


\ 

2  +  PW 


0.5 


(91) 
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which  yields  a  =  h  for  Pw  =  l.  This  equality  was  supported  by  the  depth 

measurements  in  the  lower  swash  zone  by  Kobayashi  et  al.  (1998)  who 
assumed  Pw=  1  in  Equation  89. 

The  cross-shore  velocity  U  depends  on  x  and  t  and  is  related  to  the  depth  h 
in  the  swash  zone.  The  following  relationship  between  U  and  h  may  be 
assumed  to  express  £7 as  a  function  of  h: 

U  =  ajgh+Us  (92) 

where  a  is  a  positive  constant;  and  Us  =  steady  velocity  which  is  allowed  to 
vary  with  x.  The  steady  velocity  Us  is  intended  to  account  for  offshore 

return  flow  on  the  seaward  slope  and  the  downward  velocity  increase  on  the 
landward  slope.  Holland  et  al.  (1991)  measured  the  bore  speed  and  flow 
depth  on  a  barrier  island  using  video  techniques  and  obtained  a~2  where 
the  celerity  and  fluid  velocity  of  the  bore  are  assumed  to  be  approximately 
the  same.  Tega  and  Kobayashi  (1996)  computed  wave  overtopping  of  dunes 
using  the  nonlinear  shallow-water  wave  equations  and  showed  a  ~  2  for  the 
computed  U  and  h.  As  a  result,  use  may  be  made  of  a  =  2  as  a  first  approxi¬ 
mation.  Equation  92  implies  that  the  cross-shore  velocity  U increases 
monotonically  with  the  increase  of  h  at  given  x.  Equation  92  yields  U  =  Us 
when  h  =  o,  which  may  be  acceptable  in  view  of  the  very  small  depth  in  the 

wet  and  dry  zone.  Using  Equations  89  and  92,  the  mean  (U )  and  standard 
deviation  ( )  of  the  cross-shore  velocity  U  can  be  expressed  as: 


u  =  ^{Pw9h)°'5+PwUs;  (93) 

ol  =  a2gh  -2{U-U,)(U- PJJ, )  +  Pw (0 - U,  f  (94) 

Equation  92  is  substituted  into  Equations  87  and  88  which  are  averaged 
for  the  wet  duration  using  Equation  89.  The  continuity  Equation  87  yields: 


+  U,h  =  q, 


(95) 


After  lengthy  algebra,  the  cross-shore  momentum  Equation  88  is 
expressed  as: 
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d_ 

dx 


B 


gh2  ,  g2 


=  -gSbIh-^a2ghGb(_r,) 


(96) 


with 


B  = 


2_  — 
16 


a2  + 1 ;  r  = 


3\fn  Ush 
4  g0-^ 


(97) 


The  function  Gb  (rs )  in  Equation  96  with  r  =  rs  for  simplicity  is  given  by: 


Gb  (r)  =  1  +  4nr  +  r2  for  r  >  0 


(98) 


Gft(r)  =  2exp(-r2)-r2 -l  +  Vjir[2e7/(r)  +  l]  for  r< 0  (99) 


where  erf  is  the  error  function.  The  function  Gb  increases  monotonically 
with  the  increase  of  r  and  Gb  =  o  and  1  for  r  =  -0.94  and  0.0,  respectively, 

as  shown  in  Figure  4.  For  r  <  -1.5,  Gb  —  —(l  +  \fii r  +  r2^ .  Gb  is  a  modifier  to 

the  time-averaged  shear  stress,  accounting  for  the  probability  distribution 
of  the  near-bed  velocity  in  random  waves. 


Figure  4.  Function  Gb(r)  for  wet  and  dry  zone. 
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Equation  95  and  96  are  used  to  predict  the  cross-shore  variation  of  h  and 
Us  for  assumed  qo  where  o  ,U  and  ov  are  computed  using  Equations  91, 
93  and  94,  respectively.  It  is  necessary  to  estimate  the  wet  probability  Pw 

empirically.  To  simplify  the  integration  of  the  momentum  in  Equation  96, 
the  following  formula  is  adopted: 


P„  = 


1 +  A) 


K 


A 


va 


;  A  =  q°-0  for  x<xc  (100) 
Bgh? 


where  hi  =  mean  water  depth  at  the  location  of  Pw  =  1;  n  -  empirical 
parameter  for  PW;A  =  parameter  related  to  the  wave  overtopping  and  over¬ 
flow  rate  qo  normalized  by  the  depth  hi  where  water  is  present  always.  The 
transition  from  the  wet  ( Pw  =  1)  always  zone  to  the  wet  and  dry  ( Pw  <  1) 
zone  may  be  taken  at  x  =  xSWL  where  xSWL  is  the  cross-shore  location  of  the 
still  water  shoreline  of  an  emerged  slope  (Figure  5)  or  the  seaward  edge  of  a 
submerged  crest  as  discussed  in  relation  to  Equations  85  and  86.  Equa¬ 
tion  100  is  assumed  to  be  valid  on  the  seaward  slope  and  crest  in  the  region 
of  v  <  xc  where  xc  =  landward  end  of  the  horizontal  crest  in  Figure  5. 


Figure  5.  Transition  from  wet  model  (x  <  xr )  to  wet  and  dry  model  (x  >  xswl)  for  emerged 

structure  ( Rc  >  0). 


Integration  of  Equation  96  using  Pw  given  by  Equation  100  with  h  =  hi  at 
x  =  x±  yields  h(x)  for  x±<x<xc: 


BJl  +  A)^ 


f/z. 

n-1 

1 

-1 

2  * 

zb  M  ~zb(xi)  +  ^-f  fbGbdx  (101) 
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where  Bn  =  B  (2  -  n)/ (n  - 1) ;  and  zb(x)  =  bottom  elevation  at  the  cross¬ 
shore  location  x.  The  mean  water  depth  h  at  given  x  is  computed  by  solving 
Equation  101  iteratively  where  the  bottom  friction  factor  fb  is  allowed  to 

vary  with  x  and  the  function  Gb  given  by  Equations  98  and  99  depends  on 
rs  defined  in  Equation  97.  The  empirical  parameter  n  is  taken  to  be  in  the 
range  of  1  <  n  <  2  so  that  Bn  >0.  The  formula  for  n  calibrated  using  the 
107  tests  of  wave  overtopping  and  overflow  on  a  dike  by  Farhadzadeh  et  al. 
(2007)  is  expressed  as  n  =  1.01  +  0.98[tanh(A)j  which  yields  bounds  of 
1.01  <  n<  1.99. 

The  wave  overtopping  and  overflow  rate  qQ  is  predicted  by  imposing 
Us  =  0  in  Equation  95  at  the  location  of  xc ,  namely 


Qo 


at  x  =  xc 


(102) 


where  hc  and  Pc  are  the  computed  mean  depth  h  and  wet  probability  Pw 
at  xc .  The  wave  overtopping  probability  Pa  may  be  related  to  the  wet 
probability  Pc  at  x  =  xc  where  both  Pa  and  Pc  are  in  the  range  of  0.0-1. o. 

An  empirical  relation  of  P0  =  [tanh(5Pc)|  is  fitted  for  the  107  tests  by 
Farhadzadeh  et  al.  (2007). 

On  the  slope  landward  of  the  crest,  the  wet  probability  Pw  is  assumed  to 
be  constant  and  equal  to  Pc : 

Pw  =  Pc  for  x  >  xc  (103) 


Substituting  Equation  103  into  Equation  96  and  integrating  the  resulting 
equation  from  xc  to  x,  the  mean  depth  h(x)  on  the  landward  slope  in  the 
region  of  x  >  xc  is  expressed  as: 
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(104) 
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where  the  bottom  elevation  zb  (x)  decreases  with  the  landward  increase  of 
x  in  the  region  of  x  >  xc .  Equation  104  is  solved  iteratively  to  compute  h 
at  a  given  x. 

For  assumed  qa ,  the  landward  marching  computation  of  h,  on ,  U  and  ou 

is  initiated  using  the  wet  model  in  Chapter  4  from  the  seaward  boundary 
x  =  o  to  the  landward  limit  located  at  x  =  xr  which  corresponds  to  the 
location  where  the  computed  h  or  a  becomes  negative  in  the  region  of  h 

less  than  1  cm  for  an  emerged  crest  as  shown  in  Figure  5.  For  a  submerged 
crest,  the  landward  limit  of  xr  is  taken  as  xc .  The  landward  marching 

computation  is  continued  using  the  wet  and  dry  model  in  this  section  from 
the  location  of  x  =  xSWL  where  h  =  hi  in  Equation  101  to  the  landward  end 

of  the  computation  domain  or  until  the  mean  depth  h  becomes  less  than 
0.001  cm.  Then,  the  rate  qo  is  computed  using  Equation  102.  This  landward 

computation  starting  from  qo  =  o  is  repeated  until  the  difference  between 
the  computed  and  assumed  values  of  qo  is  less  than  1  percent.  This 
convergency  is  normally  obtained  after  several  iterations.  The  computed 
values  of  h,o}],U  and  o(/  by  the  two  different  models  in  the  overlapping 
zone  of  xSWL  <  x  <  xr  (see  Figure  5)  are  averaged  to  smooth  the  transition 
from  the  wet  zone  to  the  wet  and  dry  zone. 

Kobayashi  et  al.  (2010)  compared  this  hydrodynamic  model  for  the 
impermeable  wet  and  dry  zone  with  the  107  tests  by  Farhadzadeh  et  al. 
(2007)  and  the  100  tests  conducted  by  van  Gent  (2002b)  who  measured 
the  water  depth  and  velocities  on  the  crest  and  landward  (inner)  slopes  of 
six  different  dikes.  The  agreement  was  mostly  within  a  factor  of  two  for  the 
wave  overtopping  rates  and  probabilities  as  well  as  the  water  depth, 
velocity,  and  discharge  on  the  crest  and  landward  slope  exceeded  by  2 
percent  of  the  incident  waves.  Kobayashi  et  al.  (2010a)  modified  Equations 
101  and  104  to  allow  the  integration  of  Equation  96  starting  from  an 
arbitrary  location  landward  of  the  still  water  shoreline.  This  modification 
is  necessary  for  a  berm  that  is  slanted  downward  toward  the  toe  of  a  dune. 
The  wet  probability  Pw  on  the  downward  berm  slope  is  assumed  to  be  the 
same  as  that  at  the  seaward  end  of  this  downward  slope  in  the  same  way  as 
in  Equation  103  for  the  downward  dune  slope. 
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Sediment  transport 

The  sediment  transport  model  for  the  wet  zone  in  Chapter  5  is  adjusted  for 
the  wet  and  dry  zone.  Normally  incident  waves  and  alongshore  uniformity 
are  assumed  here.  The  Gaussian  velocity  distribution  has  been  assumed  in 
Chapter  5,  whereas  U  in  the  wet  and  dry  zone  is  expressed  as  Equation  92 
along  with  the  exponential  distribution  of  h  given  by  Equation  89. 

The  movement  of  sediment  particles  is  assumed  to  occur  when  the 
instantaneous  bottom  shear  stress  given  by  0.5p  fbU2  exceeds  the  critical 
shear  stress  pg(s  - 1  )d50rpc.  as  has  been  assumed  for  Equation  49.  The 
probability  Pb  of  sediment  movement  is  then  the  same  as  the  probability 

[i0  5 

2p(s-l)d50\pc/b  1  .  Using  Equations  89  and 

92,  Pb  can  be  shown  to  be  given  by: 


Ph  =  P  for  U  >  U,, 


(105) 


pb  =  pw  exP 


Pju,b-usf 


a2gh 


for  Us<Ucb 
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w  \  cb  s 


a2gh 


for  —Us>  Ucb(107) 


where  the  upper  limit  of  Pb  is  the  wet  probability  Pw  because  no  sediment 
movement  occurs  during  the  dry  duration. 

Sediment  suspension  is  assumed  to  occur  when  the  instantaneous 
turbulent  velocity  estimated  as  (  fb  /  2)1/3  \U\  exceeds  the  sediment  fall 
velocity  wf  as  has  been  assumed  for  Equation  50.  The  probability  Ps  of 
sediment  suspension  is  then  the  same  as  the  probability  of  \U\  >  Ucs  where 
Ucs  =  wf  (2  /  fbf/ 3  ■  The  probability  Ps  is  then  given  by: 

PS=PW  for  Us  >  Ucs  (108) 


Ps  =  pw  exp 
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P  =  P  1 1  -  exp 

P  (ua  +uj 

+  exp 

pw{ua  +t02|] 

s  i 

agh 

agh  Jj 

(110) 


which  reduces  to  Equations  105-107  if  Ucs  is  replaced  by  Ucb .  If  P'>P» 
use  is  made  of  Ps  =  Pb  because  sediment  suspension  occurs  only  when 
sediment  movement  occurs. 

The  suspended  sediment  volume  Vs  per  unit  horizontal  bottom  area  is 
estimated  using  Equation  51  with  Sby  =  0  for  alongshore  uniformity.  In 
the  wet  and  dry  zone,  Vs  is  assumed  to  be  given  by: 

Vs=P,Vv(l  +  Slp  (111) 


where  Vbj=  potential  suspended  sediment  volume  on  a  horizontal  bottom 
when  Ps= 1 .  The  value  of  Vbj  is  assumed  to  be  constant  and  chosen  so  that 
the  suspended  sediment  volume  Vs  is  continuous  at  x  =  xswl  at  the 
seaward  end  of  the  wet  and  dry  zone.  The  constant  Vb/  is  based  on  the 
assumption  that  suspended  sediment  in  the  swash  zone  tends  to  remain 
suspended.  It  is  noted  that  Ps  given  by  Equations  108  -  110  decreases 
landward  with  the  decrease  of  Pw. 


Kobayashi  et  al.  (2010)  estimated  the  cross-shore  suspended  sediment 
transport  rate  qsx  is  estimated  using  Equation  52: 


=  aJJV,  ;ax 


a  +  (Sfex/tancp)°'5 


>  a 


(112) 


where  U  is  given  by  Equation  93.  The  parameter  ax  had  to  be  taken  as 
unity  in  the  zone  of  U  >  0  over  the  dune  crest  to  predict  minor  wave 
overwash.  However,  Equation  112  was  found  to  underpredict  major  wave 
overwash  in  the  three  small-scale  tests  conducted  by  Figlus  et  al.  (2009)  to 
investigate  the  transition  from  minor  to  major  wave  overwash  of  dunes 
constructed  of  fine  sand.  For  these  tests,  suspended  load  was  computed  to 
be  dominant.  To  account  for  the  wave  overtopping  rate  q0  explicitly, 
Equation  112  is  modified  as: 

q„=(a<U  +  a0U0)V,-,Uc  =  qo/h 


(113) 
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where  a0  =  empirical  parameter  with  a0= o  in  Equation  112;  and  U0  = 
onshore  current  due  to  q0 ,  which  is  significant  only  in  the  zone  of  the  very 
small  depth  h .  The  parameter  ax  is  the  same  as  in  Equation  112  without 
any  adjustment  in  the  zone  of  U  >  0 .  The  calibrated  value  for  the  three 
tests  by  Figlus  et  al.  (2009)  was  in  the  range  of  a0  =  1.3  -  1.8.  However,  the 
range  of  a0  =  0.1  -  0.5  was  necessary  for  the  minor  wave  overwash  data 
used  by  Kobayashi  et  al.  (2010a)  to  calibrate  Equations  111  and  112.  The 
accurate  prediction  of  wave  overtopping  and  overwash  is  very  difficult 
because  of  the  small  water  depth  and  large  velocity  in  the  zone  which  is 
wet  intermittently. 

The  cross-shore  bedload  transport  rate  qbx  is  estimated  using  Equation  56 
for  the  case  of  normally  incident  waves  (sin  0  =  o)  and  no  longshore  current 
[v  =  oj  where  oT  =  ou  for  sin  0  =  o  in  Equation  32.  For  this  case,  qbx  is 

given  by: 


<114) 

where  the  bottom  slope  function  Gs  (Sbx)  is  given  by  Equations  58  and  59, 
and  the  standard  deviation  o^  is  given  by  Equation  94  for  the  wet  and  dry 
zone.  The  parameter  b  in  the  wet  and  dry  zone  is  chosen  so  that  the  value 
of  qb  is  continuous  at  x  =  xswl. 

The  cross-shore  sediment  transport  rates  qsx  and  qbx  computed  for  the  wet 

zone  and  the  wet  and  dry  zone  are  averaged  in  the  overlapping  zone  of 
xswl  <  x  <  xr  f°r  the  smooth  transition  between  the  two  zones  in  the  same 

way  as  the  smooth  transition  of  h,on,U  and  o(/  as  explained  at  the  end  of 

the  “Water  depth  and  velocity”  section  in  this  chapter.  The  linear  extrapola¬ 
tion  for  the  case  of  no  overwash  given  by  Equation  60  for  scarping  is  not 
applied  now  that  the  sediment  transport  in  the  wet  and  dry  zone  is  pre¬ 
dicted.  The  continuity  equation  of  bottom  sediment  given  by  Equation  61 
with  qy  =  o  is  solved  numerically  to  obtain  the  bottom  elevation  at  the  next 

time  level. 
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9  Model  for  Permeable  Wet  and  Dry  Zone 

The  model  in  Chapter  8  is  extended  to  a  permeable  wet  and  dry  zone.  The 
extended  model  is  calibrated  and  verified  using  available  data  for  stone 
structures. 

A  number  of  time-dependent  hydrodynamic  models  for  rubble  mound 
structures  have  already  been  developed  as  reviewed  by  Losada  et  al.  (2008). 
These  numerical  models  try  to  predict  the  temporal  and  spatial  variations  of 
wave  dynamics  as  accurately  as  possible.  The  computation  time  normally 
increases  with  the  increase  of  the  resolution  and  accuracy.  The  computa¬ 
tionally  advanced  models  are  used  to  predict  hydrodynamic  variables  for 
relatively  short  durations.  To  reduce  computation  time  considerably, 
Kobayashi  et  al.  (2007b)  proposed  the  probabilistic  model  CSHORE.  The 
time-varying  wave  variables  are  expressed  using  a  probability  distribution. 
The  spatial  variations  of  the  mean  and  standard  deviation  are  computed 
using  the  time-averaged  governing  equations.  The  probabilistic  time- 
averaged  model  requires  additional  assumptions  but  its  computational 
efficiency  allows  the  calibration  of  the  model  parameters  using  a  large 
number  of  tests.  This  probabilistic  model  for  the  wet  zone  on  the  permeable 
armor  layer  is  extended  in  this  section  to  the  wet  and  dry  zone  to  predict  the 
wave  motion  above  the  still  water  level  (SWL).  The  extended  model 
provides  the  hydrodynamic  input  to  a  damage  progression  model  that 
predicts  the  slow  evolution  of  the  armor  layer  profile. 

The  movement  of  individual  stone  units  on  the  armor  layer  may  be 
computed  using  the  equation  of  motion  for  each  armor  unit  (Kobayashi 
and  Otta  1987).  The  profile  evolution  of  the  armor  layer  may  then  be 
predicted  by  computing  the  displacements  of  all  the  armor  units  (Norton 
and  Holmes  1992).  However,  this  approach  has  never  been  adopted  for 
practical  applications  probably  because  of  its  computation  time.  The 
sediment  transport  model  in  Chapter  8  is  modified  in  this  section  to 
predict  the  profile  evolution  of  the  armor  layer  in  the  same  manner  as  the 
prediction  of  the  beach  profile  evolution.  This  simple  approach  neglects 
the  discrete  nature  of  armor  stone  units  but  is  very  convenient  for  the 
prediction  of  the  armor  layer  profile  evolution  averaged  alongshore  where 
the  alongshore  averaging  reduces  the  discrete  nature. 
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Water  depth  and  velocity 

Figure  6  depicts  the  permeable  stone  layer  analyzed  in  this  section. 
Alongshore  uniformity  and  normally  incident  waves  are  assumed.  The 
cross-shore  coordinate  x  is  positive  onshore  with  x=o  at  the  offshore 
location  of  the  wave  measurement.  The  vertical  coordinate  z  is  positive 
upward  with  z=o  at  the  datum.  The  still  water  level  (SWL)  above  the 
datum  is  allowed  to  vary  in  time  during  a  storm  or  an  experiment.  The 
upper  and  lower  boundaries  of  the  permeable  stone  layer  are  located  at  z  = 
z/,(x)  and  zp(x),  respectively,  where  the  lower  boundary  is  assumed  to  be 
impermeable  to  simplify  the  analysis.  The  crest  height  Rc  is  taken 
conventionally  as  the  structure  height  above  SWL.  The  crest  location  xc  is 
defined  here  as  the  highest  and  most  landward  location.  The  wave 
overtopping  rate  is  denoted  as  q0.  The  SWL  shoreline  on  the  seaward  slope 
is  located  at  xswl.  The  mean  water  level  (MWL)  is  located  at  z  =  (s  +  rjj 

where  r|  is  the  wave  setup  above  SWL.  The  mean  water  depth  h  above  z  = 
Zb  is  given  by  h  =  (s  +  r|  -  zb  j .  The  cross-shore  location  xr  is  the  landward 
limit  of  the  time-averaged  model  in  the  wet  zone  presented  in  Chapter  6. 


Figure  6.  Transition  from  wet  model  (x<xr)  to  wet  and  dry  model  (x>xswl)  on  permeable  stone  layer. 


The  time-averaged  model  for  the  permeable  slope  developed  by  Kobayashi 
et  al.  (2007b)  has  been  modified  using  linear  wave  and  current  theory 
where  wave  overtopping  induces  onshore  current.  The  time-averaged 
continuity,  momentum,  and  wave  action  equations  are  used  to  predict  the 
cross-shore  variations  of  the  mean  U  of  the  depth-averaged  cross-shore 
velocity  U,  the  mean  q  of  the  free  surface  elevation  r|  above  SWL,  and  the 
free  surface  standard  deviation  or]  .  The  overbar  denotes  time  averaging. 

The  root-mean-square  (rms)  wave  height  is  defined  as  Hrms  =  >/8 a  . 
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Linear  progressive  wave  theory  is  used  locally  to  express  the  velocity 
standard  deviation  av  in  terms  of  ar] .  The  probability  distributions  of  r| 

and  U  are  assumed  to  be  Gaussian.  The  equivalency  of  the  time  averaging 
and  probabilistic  averaging  is  assumed  to  express  the  time-averaged  terms 
in  the  governing  equations  in  terms  of  r\,o^,U,au .  The  permeability  effects 

are  included  in  Chapter  6. 

The  landward-marching  computation  using  this  model  for  the  wet  zone  is 
continued  as  long  as  the  computed  h  and  ar]  are  larger  than  o.i  cm.  The 

end  location  of  the  computation  is  denoted  as  xr  in  Figure  6.  The  time- 
average  model  for  the  wet  zone  cannot  predict  wave  overtopping.  Conse¬ 
quently,  a  separate  model  for  the  wet  and  dry  zone  is  developed  and 
connected  with  the  model  for  the  wet  zone.  This  procedure  is  the  same  as 
that  used  in  Chapter  8.  The  time-averaged  cross-shore  continuity  and 
momentum  equations  derived  from  the  nonlinear  shallow-water  wave 
equations  on  the  permeable  slope  (Wurjanto  and  Kobayashi  1993)  are 
expressed  as: 


d 

dx 


(115) 


d 

dx 


hU2  +—h2 


rdz 
-gh 


dx 


^fb\U\U-uhwp 


(116) 


where  h  and  U  =  instantaneous  water  depth  and  cross-shore  velocity, 
respectively;  wp  =  vertical  seepage  velocity  which  is  taken  to  be  positive 
downward;  g  =  gravitational  acceleration;  Zb  =bottom  elevation  above  the 
datum  z=o;/&  =  bottom  friction  factor  which  is  allowed  to  vary  spatially; 
and  Ub  =  horizontal  fluid  velocity  at  z  =  z/>.  The  last  term  on  the  right  hand 
side  of  Equation  116  represents  the  time-averaged  flux  of  the  horizontal 
momentum  into  the  permeable  layer.  The  overbar  in  Equations  115  and 
116  for  the  wet  and  dry  zone  indicates  time  averaging  for  the  wet  duration 
only  because  no  water  exists  during  the  dry  duration. 

The  continuity  and  approximate  momentum  equations  for  the  flow  inside 
the  permeable  layer  are  expressed  as 


dx 


(117) 
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“,+PiKI  K=-sf 


ap  = 1000 
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;  Pi  =  1  p 


n  50 


A,50np 


(118) 


(119) 


where  qp  =  time-averaged  horizontal  volume  flux  in  the  permeable  layer; 
Up  =  time-averaged  horizontal  discharge  velocity;  ap  and  P1  = 

coefficients  associated  with  the  laminar  and  turbulent  flow  resistance  in 
Equation  65,  respectively;  np  =  porosity  of  the  permeable  layer;  Dn5o  = 
nominal  stone  diameter;  and  v  =  kinetic  viscosity  of  the  fluid.  Equation 
119  is  based  on  the  formula  developed  by  van  Gent  (1995)  and  calibrated 
by  Kobayashi  et  al.  (2007b).  The  resistance  component  associated  with  the 
oscillatory  flow  is  simply  neglected  in  Equation  118  which  is  solved 
analytically  to  obtain  the  discharge  velocity  Up  driven  by  the  horizontal 

pressure  gradient  due  to  q  =  (h  +  zb  ~sj  where  h  and  Zb  vary  with  x.  It  is 

noted  that  Equation  118  retains  only  the  leading  terms  in  the  horizontal 
momentum  equation  given  by  Wurjanto  and  Kobayashi  (1993).  Adding 
Equations  115  and  117  and  integrating  the  resulting  equation  with  respect 
to  x,  the  vertically  integrated  continuity  equation  is  obtained: 


hU  +  qp=q0  (120) 

where  the  wave  overtopping  rate  q0  is  defined  as  the  sum  of  the  volume 
fluxes  above  and  inside  the  permeable  layer  in  the  same  way  as  in 
Equation  62.  The  volume  flux  qp  is  estimated  as: 

9P=nA(T-zP)  <121> 

where  Pw  =  wet  probability  defined  as  the  ratio  between  the  wet  and  entire 
durations;  r\p  =  average  water  level  inside  the  permeable  layer;  and  zp  = 

elevation  of  the  impermeable  lower  boundary.  The  elevation  qp  and  zp  are 

relative  to  the  datum  z  =  o  in  Figure  6  and  [r\p  —  zp  j  is  the  thickness  of 

water  inside  the  permeable  layer.  The  elevation  qp  is  estimated  as: 


r\p  =  Pu,zb+(1-pw)zP  for  ZP>S 


(122) 
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%=pwzb+i1-pw)s  for  ZP<S  (123) 

The  upper  bound  of  r|p  for  Pw  =  l  is  the  upper  boundary  of  the  permeable 
layer  located  at  z  =  z/>  The  lower  bound  of  qp  for  Pw  =  o  is  the  higher  eleva¬ 
tion  of  the  lower  boundary  zp  of  the  permeable  layer  and  the  still  water  level 
S.  The  wet  probability  Pw  in  Equation  121  ensures  that  qp-  o  if  Pw-  o. 
Equations  121-  123  based  on  physical  reasoning  may  be  crude  but  are  used 
along  with  Equations  118  and  119  to  estimate  qp  for  the  known  h  and  Pw. 

The  momentum  flux  in  Equation  116  is  expressed  as: 

ubwp=c <-mPw(gh)°'5  wm  (124) 

with 

K  +fi1wm)wm=g  (125) 

where  am  =empirical  parameter;  and  wm  =  maximum  downward  seepage 

velocity  due  to  the  gravity  force,  obtained  by  solving  Equation  125  analy¬ 
tically.  The  seepage  velocity  wp  is  assumed  to  be  of  the  order  of  wm  or  less. 

/  — \0.5 

The  horizontal  velocity  Ub  at  z  =  z/>  is  assumed  to  be  of  the  order  of  [  ghj  . 

Equation  124  assumes  that  the  downward  flux  of  the  horizontal  momentum 
during  the  wet  duration  is  much  larger  than  the  upward  momentum  flux 
from  the  permeable  layer. 

The  cross-shore  variation  of  the  mean  water  depth  h  is  obtained  by  solving 
the  momentum  Equation  116  together  with  the  continuity  Equation  120. 

The  probability  density  function  /(7i )  in  the  wet  and  dry  zone  is  assumed  to 
be  exponential  and  given  by: 

P2  h 

f(Ji)  =  ^-exp—Pw=  for  h>  0  (126) 

with 

pu- JoXf(h)dh;h=  £h/(h)dh 


(127) 
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Equations  126  and  127  are  the  same  as  Equations  89  and  90  but  presented 
again  for  clarity.  The  wet  probability  Pw  equals  the  probability  of  the 
instantaneous  water  depth  h  >  o.  The  dry  probability  of  h  =  o  is  equal  to 
(1  -  Pw).  The  mean  water  depth  for  the  wet  duration  is  h  but  the  mean 
depth  for  the  entire  duration  is  equal  to  Pw  h  .  The  free  surface  elevation  r| 
above  SWL  is  given  by  q  =  (h  +  zb  -  S)  where  Zb  and  S  are  assumed  to  be 

invariant  during  the  averaging.  The  standard  deviations  of  77  and  h  are  the 
same  and  given  by: 


o 


n 


h 


\0.5 

2  +  PW 


(128) 


which  is  the  same  as  Equation  91. 

The  cross-shore  velocity  U  may  be  related  to  the  depth  h  in  the  wet  and  dry 
zone  in  the  same  way  as  in  Equation  92: 

U  =  ajgh+Us  (129) 

where  a  =  positive  constant  taken  as  a  =  2 ;  and  Us  =  steady  velocity 
which  is  allowed  to  vary  with  x.  The  steady  velocity  Us  is  included  to 
account  for  offshore  return  flow  on  the  seaward  slope  and  crest  and  the 
downward  velocity  increase  on  the  landward  slope.  Using  Equations  126 
and  129,  the  mean  U  and  standard  deviation  op  of  the  cross-shore 
velocity  U  can  be  expressed  as: 


U  =  fa(Pwghr+PwUs  (130) 

^^^gh-lfi-U^U-PJJ^  +  P^U-U,)2  (131) 

Equations  128, 130  and  131  express  o  ,U,ov  in  terms  of  h  ,  Pw,  and  Us 
which  vary  with  x. 

Equation  129  is  substituted  into  Equations  116  and  120  which  are  averaged 
for  the  wet  duration  using  Equation  126.  The  continuity  Equation  120 
yields: 
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3y/na  r 
- h 


gh 
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+  Ush=q;q  =  qo-q 


(132) 


where  q  =  volume  flux  above  the  permeable  layer.  After  lengthy  algebra, 
the  momentum  Equation  116  is  expressed  as: 
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with 


B  = 
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9  n 

16, 


a2  + 1 ;  rs 


Ush 

4  q  —  Ush 


(134) 


where  the  parameter  B  is  related  to  the  momentum  flux  term  on  the  left 
hand  side  of  Equation  116.  The  function  Gb(rs)  in  Equation  133  is  given  by 
Equations  98  and  99. 

Equations  132  and  133  are  used  to  predict  the  cross-shore  variation  of  h 
and  Us  for  assumed  q0.  It  is  necessary  to  estimate  the  wet  probability  Pw 
empirically.  To  simplify  the  integration  of  Equation  133,  the  following 
formula  is  adopted: 


P  = 


(i+A) 

k| 

-A 

kl 
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V  / 
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K  / 
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•  a -  9  .4 
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2 

(135) 


BgK 


where  h1  and  q,  =  mean  water  depth  and  volume  flux,  respectively,  at  the 

location  of  x  =  x2  where  Pw-i;n-  empirical  parameter  for  Pw  ;  and  A  and 
Ax  =  dimensionless  variables  related  to  q  and  qh  respectively.  The 
transition  from  the  wet  ( Pw  =  1  always)  zone  to  the  wet  and  dry  ( Pw  <1) 
zone  may  be  taken  at  Xx  =  xswl  where  xswl  is  the  cross-shore  location  of  the 
still  water  shoreline  of  an  emerged  crest  as  shown  in  Figure  6.  Equation 
135  is  assumed  to  be  valid  on  the  upward  slope  and  horizontal  crest  in  the 
region  of  x1<x<xc  where  xc  is  the  highest  and  most  landward  location  of 

the  structure. 
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Integration  of  Equation  133  for  Pw  given  by  Equation  135  starting  from 
h  =  h1  at  x  =  x1  yields  h(x) : 


Bn(  1  +  4A 


-l 


Z  (x)-z  (xj  +  f 

xl 


—aG, 


P  w 
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W° 


dx  (136) 


where  Bn  =  5(2  -  n)/(n- 1);  and  Zb(x )  =  bottom  elevation  at  the  cross-shore 
location  x.  The  mean  water  depth  h  at  given  x  is  computed  by  solving 
Equation  136  iteratively.  The  empirical  parameter  n  is  taken  to  be  in  the 
range  of  1  <  n  <  2  so  that  Bn>  o.  The  formula  for  n  for  the  impermeable 
wet  and  dry  zone  in  Chapter  8  is  adopted  and  expressed  as: 

n  =  1.01  +  0.98[tanh(^4o)]  where  1.01<n  <1.99and  Aa  —  q20  / [Bgh^ . 


On  the  downward  slope  in  the  region  of  x  >  xc,  the  wet  probability  Pw  is 
assumed  to  be  given  by: 


(137) 


where  Pc  and  qc  are  the  computed  wet  probability  Pw  and  volume  flux  q  at 
x  =  xc.  Substituting  Equation  137  into  Equation  133  and  integrating  the 
resulting  equation  from  xc  to  x,  the  mean  depth  h(x)  is  expressed  as: 
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(138) 


where  hc  is  the  computed  mean  depth  at  x  =  xc. 

The  wave  overtopping  rate  q0  is  predicted  by  imposing  Us  -  o  in  Equation 
132  at  the  crest  location  xc  such  that: 

+  qp  at  x  =  xc  (139) 


The  wave  overtopping  probability  P0  may  be  related  to  the  wet  probability 
Pc  at  x  =  xc  where  both  PQ  and  Pc  are  in  the  range  of  0.0  -  1.0.  The 
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empirical  relation  of  P0  =  [tanh(5Pc)]  for  the  impermeable  wet  and  dry 
zone  in  Chapter  8  is  adopted  to  estimate  P0. 

For  assumed  q0,  the  landward  marching  computation  of  h,on,U  and  au  is 

initiated  using  the  wet  model  in  Chapter  6  from  the  seaward  boundary  x-o 
to  the  landward  limit  located  at  x  =  xr.  The  landward  marching  computation 
is  continued  using  the  wet  and  dry  model  in  this  section  from  the  location  of 
x  =  xswl  where  h  =  h1  to  the  landward  end  of  the  computation  domain  or 

until  the  mean  depth  h  becomes  less  than  o.ooi  cm.  The  rate  q0  is 
computed  using  Equation  139  together  with  the  overtopping  probability  P0. 
This  landward  computation  starting  from  q0  =  o  is  repeated  until  the 
difference  between  the  computed  and  assumed  values  of  q„  is  less  than  1%. 
This  convergence  is  normally  obtained  after  several  iterations.  The 
computed  values  of  h,o  ,U  and  o{;  by  the  two  different  models  in  the 

overlapping  zone  of  xswl  <x<xr  (see  Figure  6)  are  averaged  to  smooth  the 
transition  from  the  wet  zone  to  the  wet  and  dry  zone. 


Farhadzadeh  et  al.  (2009)  compared  the  numerical  model  with  test 
explained  by  Kobayashi  and  de  los  Santos  (2007)  and  van  Gent  (2002a). 
The  empirical  formula  developed  using  the  tests  was  expressed  as: 


a  =a 


zb-zP 
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n  50 
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(140) 


where  the  constant  a  is  the  same  as  a  =  2  in  Equation  129  and 

[zb  -  zp )/  Dn50  is  the  local  thickness  of  the  permeable  layer  normalized  by 

the  nominal  stone  diameter. 


For  the  probability  density  function/ (h)  given  by  Equation  126,  the  water 
depth  he  corresponding  to  the  exceedance  probability  e  is  given  by: 


for  Pw  >  e 


(141) 


Using  Equation  129,  the  water  velocity  Ue  and  discharge  qe  corresponding 
to  the  exceedance  probability  e  are  expressed  as: 
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ue  =  a  Jgh~e  +us;qe  =  heUe  (142) 

The  probability  e  of  h  >he  at  given  x  is  not  directly  related  to  the 
probability  based  on  individual  overtopping  events.  The  probability  2 
percent  used  by  van  Gent  (2002b)  is  assumed  to  correspond  to  the  range 
of  e  =  0.01  -  0.02  where  Equation  141  is  not  very  sensitive  to  e  =  0.01  - 
0.02  as  long  as  the  wet  probability  Pw  is  larger  than  about  0.1.  The 
computed  values  of  he,  Ue  and  qe  in  CSHORE  are  based  one  =  0.01  where 
use  is  made  of  e  =  Pw/  1.1  if  Pw  <0.011  so  that  (Pw/e)  >1.1  in  Equation 

141.  Farhadzadeh  et  al.  (2009)  compared  the  measured  and  computed 
values  of  h2%,  U2%,  and  q2%.  The  agreement  was  mostly  within  the  factor  of 
two  but  the  hydrodynamic  variables  in  the  wet  and  dry  zone  are  difficult  to 
predict  accurately  due  to  the  small  water  depth  and  larger  velocity  during 
intermittent  wave  overtopping. 

Stone  movement 

The  sediment  transport  model  for  the  impermeable  sand  beach  in 
Chapters  5  and  8  is  modified  to  predict  the  movement  of  stone  armor  units 
on  a  coastal  structure.  The  probability  Pb  of  stone  movement  under  the 
Gaussian  velocity  U  in  the  wet  zone  is  estimated  assuming  that  the  stone 
movement  occurs  when  the  absolute  value  of  the  instantaneous  velocity  U 
exceeds  the  critical  velocity  Ucb  estimated  as: 

(7tt  =  [jVt9(S-l)D,i50]°5  (143) 

where  s  and  Dn5o  =  specific  gravity  and  nominal  diameter  of  the  stone, 
respectively;  and  Nc  =  empirical  parameter.  If  the  wave  height  Hc 
corresponding  to  Ucb  is  given  by  Hc  =  U2cb/ g ,  then  Equation  143  becomes 

Nc  =  H c  /  [(s  - 1) Dn50  ]  and  Nc  may  be  regarded  as  the  critical  stability 

number  for  the  stone  which  is  of  the  order  of  unity  (Kobayashi  et  al. 

2003).  Equation  49  is  based  on  the  critical  Shields  parameter  Wc  =  0.05 
for  the  initiation  of  sand  movement.  The  two  parameters  are  related  by 
Nc  =  2xVJfh  and  Equation  49  for  the  probability  Pi,  is  applicable  using 

>Ft  =  0.5  fbNc . 

Equation  143  is  adopted  here  and  Nc  is  calibrated  as  Nc  =  0.7  using  the 
damage  progression  tests  of  a  stone  structure  with  s  =  2.66  and  Dn50  = 

3.64  cm  conducted  by  Melby  and  Kobayashi  (1998).  The  probability  of 
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stone  suspension  is  estimated  using  Equation  50  where  the  stone  fall 
velocity  w/ is  estimated  using  wf  =  1.8  g(s-ljDn50  for  a  sphere  (e.g., 

Jimenez  and  Madsen  2003).  For  the  stone  with  s  =  2.66  and  Dn50  =  3.64 
cm,  Wf  =  1.4  m/s  and  the  computed  probability  of  suspension  of  this  stone 
is  essentially  zero.  The  stone  armor  units  are  assumed  to  move  like 
bedload  particles,  although  CSHORE  also  computes  the  suspended  stone 
transport  rate  using  the  formulas  developed  for  sand. 

The  probability  Pb  of  stone  movement  in  the  wet  and  dry  zone  is  obtained 
for  the  probability  distribution  of  U based  on  Equations  126  and  129.  The 
probability  Pb  of  stone  movement  is  assumed  to  be  the  same  as  the 
probability  of  \U\  >  Ucb  with  Ucb  given  by  Equation  143.  Then,  Pb  is  given  by 

Equations  105  - 107. 

The  time-averaged  volumetric  rate  qb  of  stone  transport  is  estimated  using 
the  formula  for  bedload  given  by  Equation  114  which  is  modified  as: 


Qbx  =  bPiPA-ofj  /[g{s- 1)1 ;  Br  = 


D 


n  50 


<1 


(144) 


where  b  =  bedload  parameter  specified  as  b  =  0.002  as  discussed  below 
Equation  57;  Gs  =  function  of  the  bottom  slope  given  by  Equations  58  and 
59;  Br  =  reduction  factor  due  to  limited  stone  availability;  m  =  empirical 
parameter;  and  ov  -  velocity  standard  deviation  representing  the  wave 

action  on  the  stone.  The  rate  qbx  becomes  negative  (offshore)  on  the  steep 
slope  with  Gs  <  o.  The  reduction  factor  Br  is  added  in  CSHORE  to  account 
for  the  thickness  (zb  -  zp )  of  the  stone  layer  where  Br=  1  if  (zb  -  zp )  >  Dn50 

and  Br  =  o  in  the  zone  of  zb  =  zp  and  no  stone.  The  computed  profile 

changes  are  found  to  be  insensitive  to  the  parameter  m  in  the  range  of  0.5  to 
2.0.  The  value  of  m  =  1.0  is  specified  in  CSHORE.  The  rate  qbx  of  stone 
transport  in  the  wet  and  dry  zone  is  also  estimated  using  Equation  144 
where  the  parameter  b  is  chosen  so  that  the  values  of  qbx  computed  for  the 
two  different  zones  are  the  same  at  the  still  water  shoreline  located  atx  = 
xswl  ■  The  computed  cross-shore  variations  of  qbx  in  the  two  zones  are 
averaged  in  the  overlapping  zone  of  xSWL  <  x  <  xr  for  the  smooth  transition 

between  the  two  zones.  The  temporal  change  of  the  bottom  elevation  zb  is 
computed  using  the  conservation  equation  of  stone  volume  in  the  same  way 
as  in  Chapter  8. 
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Farhadzadeh  et  al.  (2009)  compared  the  numerical  model  with  the  three 
damage  progression  tests  by  Melby  and  Kobayashi  (1998).  The  armor  stone 
was  placed  in  a  traditional  two-layer  thickness  with  the  seaward  slope  of 
1/2.  The  armor  stone  was  characterized  by  Dn50  =  3.64  cm,  s  =  2.66  and  np  = 
0.4  where  the  maximum  seepage  velocity  was  wm  =  8 .7  cm/s  using  Equa¬ 
tion  125.  The  thickness  of  the  armor  layer  was  7.3  cm.  The  test  duration  was 
in  the  range  of  8.5  to  28.5  hrs.  The  numerical  model  overpredicted  the 
deposited  area  below  SWL  at  the  end  of  the  test  mostly  because  it  does  not 
account  for  discrete  stone  units  dislodged  and  deposited  at  a  distance 
seaward  of  the  toe  of  the  damaged  armor  layer.  The  eroded  area  above  SWL 
was  predicted  better.  The  temporal  variation  of  the  eroded  area  Ae  was 
compared  using  damage  Se  defined  as  Se  =  Ae/  D?t50 .  The  numerical  model 

predicted  the  damage  progression  well  partly  because  the  critical  stability 
number  Nc  introduced  in  Equation  143  was  calibrated  to  be  Nc  =  0.7  for  the 
three  damage  progression  tests.  The  temporal  variations  of  Se  computed  for 
Nc  =  0.7  and  0.6  were  fairly  sensitive  to  Nc.  The  simple  criterion  of  stone 
movement  based  on  Equation  143  may  be  improved  so  as  to  predict  the 
damage  progression  more  accurately. 
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10  Sensitivity  Tests,  Calibration  and  Model 
Validation  to  Field  Data 

Introduction 

The  CSHORE  model  input  specification  includes  eight  sediment  related 
parameters  influencing  the  calculated  sediment  transport  and  profile 
evolution.  Four  of  these  parameters  pertain  to  the  physical  characteristics  of 
the  sediment  comprising  the  beach  profile  and  therefore  should  not  be 
adjusted  for  model  calibration  purposes.  These  parameters  are  the  median 
sediment  diameter  D50  (mm),  the  sediment  fall  velocity  in/ (m/s),  the 
sediment  specific  gravity  s,  and  the  sediment  maximum  slope  tan  &.  The 
free  parameters  include  the  empirical  bedload  parameter  b,  the  empirical 
suspended  load  parameter  a,  and  two  sediment  suspension  efficiency 
parameters  eB  and  e/for  sediment  suspension  due  to  energy  dissipation 
from  wave  breaking  and  bottom  friction  respectively.  Kobayashi  et  al. 
(2008a)  calibrated  the  bedload  parameter  b  =  0.002  using  20  water  tunnel 
tests  of  Ribberink  and  Al-Salem  (1994),  the  four  large-scale  flume  tests  of 
Dohmen- Janssen  and  Hanes  (2002)  and  the  24  sheet  flow  test  by  Dohmen- 
Janssen  et  al.  2002).  Based  on  three  small-scale  equilibrium  profile  tests 
conducted  by  Kobayashi  et  al.  (2005),  the  suspended  load  parameter  a  was 
estimated  on  the  order  of  0.2.  Typical  values  for  the  suspension  efficiency 
parameters  eB  =  0.005  and  e/=o.oi  have  be  assigned  although  these  values 
are  uncertain  (Kobayashi  2007a). 

In  this  chapter,  the  CSHORE  model  is  applied  to  a  suite  of  52  cases  of 
storm-induced  beach  profile  change  measured  in  the  field.  These  data  sets 
encompass  seven  sites  along  the  Atlantic  coast  and  two  sites  along  the 
Pacific  coast  and  involve  several  storm  events.  An  initial  suite  of  sensitivity 
simulations  was  performed  to  quantify  the  relative  influence  and  role  of 
each  of  the  model  parameters  on  the  predicted  profile  change.  CSHORE 
was  calibrated  for  three  east  and  two  west  coast  cases  with  y  set  at  the 
default  value  of  0.5.  Based  on  these  calibrations,  new  default  sediment 
parameter  values  were  estimated.  The  full  suite  of  field  cases  was 
simulated  using  the  default  parameter  values. 
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Sensitivity  tests 

To  determine  the  sediment  parameter  sensitivity,  CSHORE  was  run  with 
the  initial  default  values  for  all  but  one  of  the  empirical  inputs.  The  variable 
was  systematically  varied  over  a  reasonable  range  of  values.  The  other 
parameter  values  were  held  constant  at  their  respective  default  values. 

Bedload  parameter  b 

Figures  7  through  11  illustrate  the  sensitivity  of  CSHORE  results  to  the 
specified  value  of  the  bedload  parameter  b.  The  six  cases  illustrated  are  a 
high  quality  and  representative  data  subset  used  in  the  sensitivity  analysis. 
The  bedload  parameter  controls  the  magnitude  of  the  computed  bedload 
sediment  transport  which  is  prescribed  in  CSHORE  as  onshore-directed. 
Therefore,  as  seen  in  the  figures,  increasing  the  value  of  b  results  in  a 
reduction  of  erosion  in  the  swash  zone  and  for  larger  values,  a  buildup  of 
sediment  in  the  foreshore  is  observed.  As  the  value  of  b  decreases  the  bar 
feature  is  positioned  further  offshore  and  is  less  pronounced,  again,  the 
result  of  reduced  onshore  sand  transport  counter  acting  offshore  directed 
suspended  sand  transport.  For  these  test  cases,  values  of  b  greater  than 
about  0.005  result  in  unrealistic  accretion  of  the  beach  berm  feature. 

Suspended  load  parameter  a 

Figures  12  through  16  illustrate  the  sensitivity  of  CSHORE  results  to  the 
specified  value  of  the  suspended  load  parameter  a.  The  suspended  load 
parameter  accounts  for  the  onshore-directed  wave-related  sediment 
transport  through  a  reduction  in  the  offshore-directed  sediment  flux  due  to 
undertow.  The  undertow  driven  erosion  of  the  foreshore/swash  zone  and 
volume  of  sand  deposited  in  the  bar  feature  is  unmitigated  by  the  effect  of 
wave  transport  as  the  value  of  a  approaches  one.  Model  results  are  relatively 
sensitive  to  the  specified  value  of  the  suspended  load  parameter  a. 

Suspension  efficiency  parameter  ee 

Figures  17  through  21  illustrate  the  sensitivity  of  CSHORE  results  to  the 
specified  value  of  the  suspension  efficiency  parameter  eB.  The  suspension 
model  is  predicated  on  the  importance  of  energy  dissipation  in  sediment 
entrainment,  and  the  suspension  efficiency  parameter  eB  is  the  primary  free 
parameter  controlling  magnitude  of  transport.  Increasing  the  value  of  eB 
increases  the  suspended  sediment  concentration  in  the  breaking  region  and 
increases  the  estimated  erosion  on  the  foreshore  and  deposition  in  the 
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DUO 

Dewey  Profile  140 


Figure  7.  Case  D140  -  Sensitivity  to  bedload  parameter  b. 


R236 

Rehoboth  Profile  236 


Figure  8.  Case  R236  -  Sensitivity  to  bedload  parameter  b. 
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HAL74 

Oceancity  MD  profile  74 


Figure  9.  Case  HAL74  -  Sensitivity  to  bedload  parameter  b. 


DM580 
Del  Mar,  CA 


Figure  10.  Case  DM580  -  Sensitivity  to  bedload  parameter  b. 
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OS 1000 
Oceanside,  CA 


Figure  11.  Case  OSIOOO  -  Sensitivity  to  bedload  parameter  b. 
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Figure  12.  Case  D140  -  Sensitivity  to  suspended  load  parameter  a. 


ERDC/CHL  TR-12-22 


59 


R236 

Rehoboth  Profile  236 


Figure  13.  Case  R236  -  Sensitivity  to  suspended  load  parameter  a. 
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Oceancity  MD  profile  74 


Figure  14.  Case  HAL74  -  Sensitivity  to  suspended  load  parameter  a. 
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DM580 
Del  Mar,  CA 


Figure  15.  Case  DM580  -  Sensitivity  to  suspended  load  parameter  a. 


OS 1000 
Oceanside,  CA 


Figure  16.  Case  0S1000  -  Sensitivity  to  suspended  load  parameter  a. 
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DUO 

Dewey  Profile  140 


Figure  17.  Case  D140  -  Sensitivity  to  efficiency  parameter  eB. 
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Rehoboth  Profile  236 


Figure  18.  Case  R236  -  Sensitivity  to  efficiency  parameter  eB. 
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HAL74 

Oceancity  MD  profile  74 


Figure  19.  Case  HAL74  -  Sensitivity  to  efficiency  parameter  eB. 
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Figure  20.  Case  DM580  -  Sensitivity  to  efficiency  parameter  eB. 
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OSIOOO 
Oceanside,  CA 


Figure  21.  Case  OSIOOO  -  Sensitivity  to  efficiency  parameter  eB. 


nearshore  bar  feature.  Conversely,  a  decrease  in  the  value  of  eB  results  in 
less  foreshore  erosion  and  a  more  subtle  bar  feature.  It  appears  that  the 
CSHORE  results  may  be  more  sensitive  to  changes  in  the  value  of  eB  than 
they  are  to  changes  in  the  value  of  a.  Although  the  parameters  are 
associated  with  different  physical  mechanisms,  both  parameters  tend  to 
influence  the  model  results  in  a  similar  way. 

Suspension  efficiency  parameter  et 

As  shown  in  Kobayashi  (2008a)  the  suspended  sediment  concentration 
and  transport  is  driven  by  both  wave  breaking  dissipation  and  dissipation 
in  the  bottom  boundary  layer.  The  magnitude  of  suspended  sediment 
concentration  due  to  bottom  friction  is  controlled  through  the  efficiency 
parameter  e/.  Figures  22  through  26  illustrate  the  sensitivity  of  CSHORE 
results  to  the  specified  value  of  the  suspension  efficiency  parameter  due  to 
bottom  friction.  In  a  manner  similar  to  the  breaking  efficiency,  increasing 
the  value  of  e/ tends  to  increase  the  estimated  foreshore  erosion.  Although 
it  may  not  be  evident  in  this  effort  the  breaking  dissipation  and  bottom 
friction  dissipation  can  be  important  suspension  mechanisms  in  different 
regions  of  the  profile,  and  therefore  increases  in  these  two  values  can 
result  in  divergent  profile  evolution. 
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D140 

Dewey  Profile  140 


Figure  22.  Case  D140  -  Sensitivity  to  efficiency  parameter  et. 
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Rehoboth  Profile  236 


Figure  23.  Case  R236  -  Sensitivity  to  efficiency  parameter  et. 
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HAL74 

Oceancity  MD  profile  74 


Figure  24.  Case  HAL74  -  Sensitivity  to  efficiency  parameter  et. 
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Figure  25.  Case  DM580  -  Sensitivity  to  efficiency  parameter  et. 
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Figure  26.  Case  OSIOOO  -  Sensitivity  to  efficiency  parameter  et. 


Summary 

CSHORE  results  are  more  sensitive  to  changes  in  the  values  of  the 
suspended  load  parameter  a  and  the  suspension  efficiency  parameter  eB 
than  they  are  to  reasonable  changes  in  the  value  of  the  bedload  parameter 
b  or  the  suspension  efficiency  parameter  e/.  It  is  noted  that  the  suspended 
load  parameter  a  and  the  suspension  efficiency  parameter  eB  are 
associated  with  different  physical  processes  but  tend  to  influence  the 
model  results  in  the  same  way.  That  is,  increasing  the  value  of  either  a  or 
eB  will  act  to  increase  the  amount  of  erosion  at  the  foreshore  and  the 
volume  of  sand  in  the  nearshore  bar  feature.  The  suspension  efficiency 
parameters  eB  and  e/  influence  the  magnitude  of  the  calculated  suspended 
sediment  transport,  and  the  suspended  load  parameter  a  controls  the 
balance  between  suspended  sediment  transport  directed  offshore  due  to 
under  tow  and  onshore  due  to  the  positive  correlation  between  the  time 
varying  cross-shore  velocity  and  suspended  sediment  concentration.  The 
bedload  parameter  b  influences  the  magnitude  of  the  onshore  directed 
bedload  sediment  transport .  Increasing  the  value  of  b  will  increase  the 
onshore  bedload  transport  which  will  tend  to  decrease  the  volume  in  the 
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nearshore  bar  and  reduce  foreshore  erosion  or,  at  large  values  of  b, 
produce  foreshore  deposition. 

CSHORE  calibration 

Having  identified  the  role  and  relative  sensitivity  of  the  CSHORE  input 
parameters,  appropriate  default  parameter  values  to  serve  as  a  reasonable 
starting  point  for  project  specific  CSHORE  calibration  are  developed 
herein.  In  the  case  that  calibration  data  for  CSHORE  does  not  exist, 
predictions  maybe  made  using  default  parameter  settings.  Consequently, 
the  approach  used  in  the  present  study  was  to  perform  a  detailed 
calibration  for  a  limited  number  of  field  cases  (three  Atlantic  coast  cases 
and  two  Pacific  coast  cases)  and  default  parameter  values  were  estimated. 
The  predictive  skill  of  CSHORE  was  tested  by  running  the  model  against 
the  larger  suite  of  52  field  cases  using  the  new  default  parameter  value 
settings.  Evaluation  of  the  model  performance  focuses  on  simulation  of 
erosion  processes  at  the  foreshore  and  dune  with  less  emphasis  on  the 
submerged  portion  of  the  beach  profile. 

Overlap  in  effect  between  the  CSHORE  input  parameters  tends  to 
complicate  the  model  calibration  process;  in  that  the  best  fit  value  for  one 
parameter  depends  to  some  extent  on  the  values  assigned  to  the  other 
parameters.  As  a  result,  a  family  of  different  parameter  value  settings  is 
possible,  each  leading  to  a  satisfactory  calibration  result.  In  the  present 
calibration  effort  iterative  simulations  were  performed  between  optimized 
values  of  the  breaking  efficiency  parameter  and  the  suspended  load 
parameter.  In  these  simulations  optimization  was  based  on  minimizing  the 
root  mean  square  error  between  the  measured  and  computed  post-storm 
beach  profile  across  the  entire  profile  (computational  domain).  After 
optimum  values  for  the  breaking  efficiency  and  suspended  load  parameters 
were  determined  the  value  of  the  bedload  parameter  was  optimized.  Next 
the  value  of  the  bottom  friction  efficiency  parameter  e/was  optimized. 

Suspension  efficiency  parameter  eg 

The  value  of  suspension  efficiency  parameter  eB,  was  the  first  model 
parameter  to  be  optimized  in  the  calibration  process.  The  suspension 
efficiency  parameter  eB  optimized  at  values  ranging  from  0.004  to  0.006 
for  the  Atlantic  coast  calibration  cases.  For  the  Pacific  coast  calibration 
cases  the  suspension  efficiency  parameter  eB  optimized  at  values  ranging 
from  0.01  to  0.02.  Default  values  for  the  suspension  efficiency  parameter 
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were  set  at  0.005  for  Atlantic  coast  applications  and  twice  that,  0.01,  for 
Pacific  coast  applications.  Again,  the  differences  in  the  Atlantic  and  Pacific 
coasts  tests  are  poorly  understood,  and  the  larger  parameter  value  for 
Pacific  coast  applications  relates  to  the  requirement  for  larger  suspended 
load  sediment  transport  and  enhanced  erosion. 

Suspended  load  parameter  a 

The  suspended  load  parameter  a  optimized  at  values  ranging  from  0.1  to 
0.2  for  Atlantic  coast  applications.  The  value  of  the  suspended  load 
parameter  for  the  Pacific  coast  calibration  cases  ranged  between  0.4  and 
0.6.  The  new  default  value  for  the  suspended  load  parameter,  a,  was  set  at 
0.2  for  Atlantic  coast  applications  and  0.5  for  the  Pacific  coast.  Although 
the  difference  in  parameter  values  is  not  well  understood,  a  distinguishing 
difference  between  the  Atlantic  coast  and  Pacific  coast  calibration  cases  is 
the  relative  extent  of  upper  beach  erosion.  This  difference  stems  from  the 
comparatively  mild  pre-storm  foreshore  slopes  in  the  Pacific  coast  cases 
(approximately  1  on  35)  as  compared  to  the  pre-storm  foreshore  slopes  in 
the  Atlantic  coast  cases  (approximately  1  on  18).  Whereas,  the  post-storm 
foreshore  slopes  are  nearly  the  same  for  the  Pacific  and  Atlantic  coast 
cases  at  1  on  21  and  1  on  16,  respectively.  To  accommodate  prediction  of 
the  larger  upper  beach  erosion  observed  in  the  Pacific  coast  calibration 
cases,  the  proportion  of  offshore  directed  suspended  sediment  transport 
had  to  increase  relative  to  the  onshore  directed  bedload  sediment 
transport,  which  was  achieved  in  part  with  larger  parameter  values  for  the 
suspended  load  parameter  a. 

Bedload  parameter  b 

The  calibration  process  revealed  that  satisfactory  predictions  were 
consistently  obtained  with  a  bedload  parameter  value  between  0.001  and 
0.002.  The  new  default  for  both  the  Atlantic  and  Pacific  coast  applications 
is  set  at  the  preliminarily  assigned  default  value  of  0.002.  It  is  noted  that 
the  bedload  parameter  is  expected  to  be  instrumental  in  the  accretionary 
phase  of  profile  development.  The  short-term  erosional  field  cases 
presented  herein  are  ill-suited  to  the  calibration  of  longer-term  evolution, 
and  additional  data  may  alter  this  conclusion. 

Suspension  efficiency  parameter  ef 

The  value  of  the  suspension  bottom  friction  efficiency  e/  was  selected 
following  the  other  parameter  optimizations.  The  procedure  indicates  that 
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changes  in  the  value  assigned  to  this  model  parameter  had  very  little 
influence  on  the  predicted  beach  profile  response.  The  default  value  for 
this  model  parameter  remained  unchanged  at  o.oi  for  the  Atlantic  coast 
cases  and,  consistent  with  the  breaking  efficiency  parameter,  was  doubled 
for  the  Pacific  coast  cases  to  0.02. 

Summary 

CSHORE  model  results  for  each  of  the  six  calibration  cases  are  illustrated  in 
Figures  27  through  31.  Each  figure  plots  the  measured  initial  and  final  beach 
profile  together  with  the  CSHORE  predicted  post-storm  profile  based  on  the 
new  default  parameter  values  developed  as  part  of  the  present  effort.  The 
Atlantic  coast  cases  are  shown  in  Figures  27  through  29  whereas  the  Pacific 
coast  cases  are  shown  in  Figures  30  and  31. 


D140 

Dewey  Profile  140 


Figure  27.  Case  D140. 
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R236 

Rehoboth  Profile  236 


Figure  28.  Case  R236. 
HAL74 

Oceancity  MD  profile  74 


Figure  29.  Case  HAL74. 
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DM580 
Del  Mar,  CA 


Figure  30.  Case  DM580. 
OSIOOO 
Oceanside,  CA 


Figure  31.  Case  OSIOOO. 
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Field  case  studies 

Based  on  the  Atlantic  coast  calibration  effort  new  default  parameter  values 
were  estimated  as  follows: 

•  b  =  0.002 

•  a  =  0.2 

•  eB  =  0.005 

•  e/=  0.01 


Using  these  default  values,  CSHORE  simulations  were  performed  for  a 
suite  of  fifty  (50)  Atlantic  coast  field  cases  involving  seven  different  sites 
and  six  different  storm  sequences.  The  field  cases  selected  for  this  study 
are  the  same  data  that  were  used  to  document  the  validation  of  SBEACH 
(Wise  et  al.  1996)  and  represent  the  best  and  most  readily  available 
records  of  storm-induced  profile  response  for  application  and  evaluation 
of  CSHORE.  The  field  case  studies  include  Ocean  City,  MD  (a  series  of 
storms  from  October  1991-  January  1992);  Manasquan/Point  Pleasant 
Beach,  NJ  (March  1984  storm)  Dewey  Beach/Rehoboth  Beach,  DE 
(December  1992  storm);  and  Debidue  Beach/Myrtle  Beach,  SC  (Hurricane 
Hugo,  September  1989). 

Based  on  the  Pacific  coast  calibration  effort  new  (Pacific  coast)  default 
parameter  were  estimated  as: 

•  b  =  0.002 

•  a  =  0.5 

•  eB  =  0.01 

•  e/=  0.02 


Using  these  default  values,  CSHORE  simulations  were  performed  for  the 
three  available  field  cases  in  southern  California  (Del  Mar  and  Oceanside). 

Measuring  model  performance 

The  intent  of  this  portion  of  the  study  is  to  illustrate  the  capability  of 
CSHORE  to  reproduce  storm-induced  beach  profile  change  measured  in 
the  field  and  to  evaluate  model  performance  through  comparison  of 
measured  and  predicted  beach  profile  response  with  emphasis  on  the  dune 
and  foreshore.  The  measures  of  performance  were  selected  based  on  the 
types  of  information  that  field  engineers  require  in  the  assessment  of 
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storm  damage.  These  measures  include  volume  change  at  the  foreshore, 
recession  of  a  specified  contour,  and  landward  storm  intrusion,  see 
Figure  32  for  a  schematic. 


Figure  32.  Definition  of  model  performance  measures  at  the  foreshore. 


Volume  change.  Volume  change  at  the  foreshore  is  computed  as  the  net 
gain  or  loss  in  volume  per  unit  beach  width  above  some  specified  elevation. 
In  this  study,  the  reference  elevation  for  volume  calculations  was  set  to  the 
National  Geodetic  Vertical  Datum  (NGVD)  which  approximately  corre¬ 
sponds  to  mean  sea  level  at  the  selected  field  sites.  The  measure  of  volume 
change  is  useful  in  identifying  how  well  CSHORE  predicts  erosion  of  the 
subaerial  beach,  which  is  important  in  defining  volume  requirements  for 
beach  nourishment  operations  and  estimating  damages  to  upland  struc¬ 
tures.  However,  in  the  available  data  set,  volume  change  is  not  always  a 
reliable  indicator  of  model  accuracy  in  the  prediction  of  storm  erosion. 
Some  of  the  post-storm  surveys  occurred  well  after  the  end  of  the  erosive 
wave  conditions  and  include  berm  features  that  indicate  significant  post 
storm  beach  profile  recovery.  In  these  cases  calculation  of  volume  change 
based  on  the  pre-  and  post-storm  profiles  underestimates  the  erosion 
generated  by  the  storm  and  predicted  by  CSHORE.  Consequently,  the 
presence  of  a  post-storm  recovery  berm  must  be  considered  when 
evaluating  model  performance  using  this  measure. 

Contour  recession.  Measuring  the  horizontal  recession  of  a  specified 
contour  provides  an  assessment  of  how  well  CSHORE  predicts  the  change 
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of  a  given  morphological  feature  such  as  the  dune  or  berm.  Contour 
recession  also  provides  information  on  the  loss  of  beach  width  resulting 
from  a  storm.  For  this  study,  a  representative  beach  contour  located 
between  the  reference  elevation  (NGVD)  and  the  dune  crest  was  identified 
for  each  field  case.  Recession  of  the  selected  contour  from  its  initial 
position  on  the  pre-storm  profile  was  determined  from  the  data  and 
compared  with  simulation  results. 

Landward  storm  intrusion.  Another  measure  of  beach  profile  change 
at  the  foreshore  is  the  extent  of  landward  storm  intrusion.  Landward 
storm  intrusion  is  defined  as  the  landward  limit  of  a  specified  depth  of 
profile  erosion.  In  this  study  an  erosion  depth  of  0.3  m  was  employed  for 
all  cases.  The  measure  of  storm  intrusion  distance  is  often  used  to  quantify 
the  landward  extent  of  potential  structural  damage  caused  by  erosion  and 
undermining.  All  measures  of  landward  storm  intrusion  presented  in  this 
report  are  referenced  to  the  horizontal  position  of  the  vertical  datum 
(NGVD)  on  the  initial  beach  profile 

Ocean  City,  MD 

Between  late  October  1991  and  early  January  1992  the  beach  at  Ocean  City, 
MD  was  affected  by  a  series  of  severe  storms.  The  storms,  which  included 
the  30  October  1991  “Halloween  storm,”  11  November  1991  storm,  and  4 
January  1992  storm,  resulted  in  large-scale  erosion  of  the  constructed 
foreshore  and  dune  at  Ocean  City.  A  major  beach  nourishment  project  had 
recently  been  completed  at  Ocean  City,  the  result  of  a  joint  effort  of  the  U.S. 
Army  Corps  of  Engineers,  the  state  of  Maryland,  Worcester  County  and  the 
Town  of  Ocean  City  (Grosskoph  and  Stauble  1993).  The  data  used  in  this 
case  study  was  collected  as  part  of  the  project  monitoring  effort,  which 
included  collection  of  beach  profile  data  by  sled  surveys  and  the  measure¬ 
ment  of  local  wave  conditions  and  water  levels  at  an  offshore  gauge.  These 
monitoring  data  provided  a  comprehensive  set  of  high-quality  data  for 
application  of  CSHORE.  Figures  33-35  present  the  wave  and  water  elevation 
time  histories  recorded  by  the  offshore  gauge  and  used  as  input  to  CSHORE 
for  each  storm.  The  wave  gauge  was  located  outside  the  surf  zone  at  a  depth 
of  approximately  10  m.  Analysis  of  sediment  samples  collected  as  part  of  the 
monitoring  program  indicated  a  representative  median  sediment  grain  size 
of  0.35  mm.  Additional  details  on  the  beach  nourishment  project  and 
monitoring  program  at  Ocean  City  are  provided  by  Stauble  et  al.  (1993). 
Impacts  of  the  storms  at  Ocean  City  are  discussed  in  a  special  issue  of  Shore 
and  Beach  (Vol.  61  No.  1, 1993). 


Welter  Elevation,  m  (NGVD) 
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Time,  hr 

Figure  35.  Wave  height,  wave  period,  and  water  elevation  time-histories  for  the  4  January 

1992  storm,  Ocean  City,  MD. 

Halloween  storm.  As  shown  in  Figure  33,  the  Halloween  storm  produced 
large  waves  and  high  water  levels  at  Ocean  City  for  approximately  four  days. 
The  maximum  measured  significant  wave  height  was  about  3  m  and  the 
peak  water  level  was  measured  at  just  over  1.5  m  NGVD.  Data  at  six  beach 
profile  survey  lines  were  available  to  model  the  beach  response  to  the  storm. 
Pre-storm  profile  surveys  were  collected  in  June  1991  and  post-storm 
surveys  were  collected  shortly  after  the  storm.  The  June  surveys  were 
assumed  to  reasonably  characterize  the  condition  of  the  upper  beach  profile 
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prior  to  the  storm,  as  no  other  significant  erosion  events  occurred  between 
June  and  October.  The  storm-induced  volume  change  at  each  of  the  six 
profiles  is  shown  in  Figure  36,  where  it  is  seen  that  five  of  the  profiles 
experienced  a  net  loss  of  volume  and  one  profile  experienced  a  net  gain  in 
volume.  The  lack  of  mass  conservation  at  the  various  profiles  indicates  that 
material  was  redistributed  by  longshore  processes  between  the  surveys.  It  is 
believed  that  most  of  the  sediment  redistribution  was  confined  to  the 
subaqueous  portion  of  the  profile  and  that  changes  on  the  upper  beach 
profile  and  dune  were  dominated  by  cross-shore  sediment  transport 
processes  during  the  storm.  Consequently,  model  results  are  expected  to 
produce  better  agreement  on  the  upper  beach  profile  than  on  the  offshore 
portion  of  the  profile. 


Figure  36.  Conservation  of  mass  between  beach  profile  surveys  for  the 
Halloween  storm,  Ocean  City,  MD. 


CSHORE  results  for  the  Halloween  storm  are  provided  in  Figures  3 7 
through  42.  The  beach  profile  response  to  this  storm  is  reasonably  well 
predicted  by  CSHORE  simulations  for  five  of  the  six  profiles.  For  profiles  3 7, 
63,  and  74  predicted  erosion  above  datum  compares  favorably  with  the 
measurements.  For  profiles  45  and  56,  predicted  erosion  above  the  datum  is 
slightly  less  than  indicated  in  the  measurements.  For  profile  103  the 
predicted  erosion  largely  exceeds  the  erosion  indicated  in  the  measure¬ 
ments.  The  measured  nearshore  recovery  bar  or  beach  step  present  in  the 
post-storm  profiles  is  not  indicated  in  the  CSHORE  results. 
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Figure  37.  CSHORE  simulation  of  the  Halloween  storm  for  Profile  37,  Ocean  City,  MD. 


HAL45 

Oceancity  MD  profile  45 


Figure  38.  CSHORE  simulation  of  the  Halloween  storm  for  Profile  45,  Ocean  City,  MD. 
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Oceancity  MD  profile  56 


ire  39.  CSHORE  simulation  of  the  Halloween  storm  for  Profile  56,  Ocean  City,  MD. 


HAL63 

Oceancity  MD  profile  63 


Figure  40.  CSHORE  simulation  of  the  Halloween  storm  for  Profile  63,  Ocean  City,  MD. 
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Figure  41.  CSHORE  simulation  of  the  Halloween  storm  for  Profile  74,  Ocean  City,  MD. 
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Figure  42.  CSHORE  simulation  of  the  Halloween  storm  for  Profile  103,  Ocean  City,  MD. 
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Figures  43  -  45  provide  plots  of  the  quantitative  model  performance 
measures  described  above.  Figure  71  shows  that  although  the  calculated 
volumetric  erosion  above  NGVD  is  generally  under  estimated,  the 
measured  profile  to  profile  trends  were  followed  except  for  profile  103. 

Figure  44  shows  the  measured  and  calculated  recession  of  the  1.5  m 
contour.  Here  it  is  seen  that  for  four  of  the  six  profiles  the  predicted 
recession  of  the  1.5  m  contour  is  within  about  five  or  six  meters  of  the 
measured  with  a  bias  for  under  prediction.  The  comparisons  for  profiles 
56  and  103  are  not  as  good  with  an  under  prediction  of  about  12  m  at 
profile  56  and  an  over  prediction  of  about  17  m  at  profile  103.  It  is  noted 
that  profile  103  gained  approximately  100  cu  m/m  between  the  pre-  and 
post-storm  surveys,  which  indicates  that  this  profile  was  heavily 
influenced  by  non-uniform  longshore  sediment  transport  processes  and 
may  explain  why  the  CSHORE  predictions  are  poorest  at  this  profile. 
Figure  45  shows  the  comparison  between  the  measured  and  calculated 
storm  intrusion  of  the  0.3-m  depth  of  erosion.  The  difference  between  the 
calculated  and  measured  storm  intrusion  distance  is  less  than  about  10  m 
except  for  at  profile  103,  which  is  considered  reasonably  good  agreement. 


Figure  43.  Comparison  of  measured  and  calculated  volume  change  above 
NGVD  for  the  Halloween  storm,  Ocean  City,  MD. 


Storm  Intrusion  (m)  era'  Contour  Recession  (m) 
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Recession  of  +1 .5  m  Contour 
Ocean  City,  Maryland  (Halloween  Storm) 


jre  44.  Comparison  of  measured  and  calculated  recession  of  the  1.5  m 
contour  for  the  Halloween  storm,  Ocean  City,  MD. 


Landward  Storm  Intrusion 
Ocean  City,  Maryland  (Halloween  Storm) 


Figure  45.  Comparison  of  measured  and  calculated  storm  intrusion  for  the 
Halloween  storm,  Ocean  City,  MD. 
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November-January  storms.  The  November-January  (NJ)  cases 
involved  a  combined  simulation  of  the  11  November  1991  storm  and  the  4 
January  1992  storm.  It  was  necessary  to  combine  the  two  storms  because 
beach  profile  data  were  not  available  between  the  two  events.  The  wave 
height  and  water  level  time  histories  are  shown  in  Figures  34  and  35.  As 
indicated  in  the  figures  the  duration  of  both  storms  was  approximately 
three  days.  The  November  storm  involved  a  sustained  maximum 
significant  wave  height  of  approximately  3  m  for  about  10  hours  and  a 
peak  total  water  elevation  of  1.2  m  NGVD.  The  January  storm  was 
somewhat  more  energetic  than  the  November  storm  with  significant  wave 
heights  exceeding  3  m  for  a  15  hour  duration  peaking  at  approximately 
4  m.  The  peak  total  water  elevation  was  measured  at  2  m.  The  January 
storm  is  the  most  energetic  of  the  three  Ocean  City  storms  included  in  this 
study  in  terms  of  maximum  significant  wave  height  and  peak  total  water 
elevation  and  this  storm  also  produced  the  most  damage  in  terms  of 
erosion  on  the  upper  beach  profiles. 

For  the  NJ  simulations  the  beach  profiles  measured  after  the  Halloween 
storm  served  as  the  pre-storm  profiles  and  beach  profile  surveys  collected 
soon  after  the  January  storm  provided  the  post-storm  profiles  used  for 
comparison  with  the  model  results.  Figure  46  illustrates  the  degree  of 
mass  conservation  between  the  pre-  and  post-storm  surveys.  Here  it  is 
seen  that  profiles  37  and  103  gained  considerable  mass  between  the  survey 
interval  whereas  mass  is  nearly  conserved  at  the  other  profile  lines. 

The  measured  and  CSHORE  predicted  beach  profile  response  for  the  NJ 
storm  series  is  shown  in  Figures  47  through  52.  These  figures  indicated  that 
the  NJ  storm  series  produced  significant  dune  overwash,  a  reduction  in 
dune  crest  elevation  and  a  landward  translation  of  the  dune  crest  at  profiles 
45,  63,  74  and  103.  Although  the  CSHORE  simulations  for  these  profiles 
included  hydrodynamic  overtopping  of  the  dune  crest  the  predicted  dune 
crest  lowering  is  far  less  than  indicated  in  the  measurements  and  no 
landward  translation  of  the  dune  crest  was  predicted.  Sensitivity  testing 
showed  that  the  predicted  dune  lowering  could  be  increased  by  increasing 
the  value  of  the  breaker  ratio  y.  However,  use  of  a  larger  breaker  ratio 
specification  in  other  simulations  resulted  in  predicted  dune  crest  lowering 
in  cases  where  the  measurements  indicate  no  dune  crest  lowering.  These 
findings  indicate  that  improvements  are  needed  in  CSHORE  in  the  area  of 
dune  overtopping  and  overwash  processes. 


Elevation  [m]  NGVD  I  (§.:  |  Elevation  [m)  NGVD 
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NJ103 

Oceancity  MD  profile  103  Nov  thru  Jan 


Figure  52.  CSHORE  simulation  of  the  NJ  storm  series  for  Profile  103,  Ocean  City,  MD. 


The  quantitative  model  performance  measures  are  illustrated  in  Figures  53 
through  55.  In  Figure  53  it  is  seen  that  the  measured  and  predicted  volume 
change  above  NGVD  compare  favorably  at  profiles  45,  63,  and  103  even 
though  the  predicted  and  measured  profiles  do  not.  At  profile  56  the 
predicted  volume  change  above  NGVD  as  well  as  the  predicted  profile 
compare  well  with  the  measurements.  At  profile  37  the  measurements 
indicate  a  net  gain  of  approximately  8  cu  m/m  whereas  the  CSHORE 
simulation  estimated  a  net  loss  of  more  than  20  cu  m/m.  At  profile  74 
CSHORE  underestimates  volume  loss  above  NGVD  by  nearly  25  cu  m/m 
presumably  a  result  of  the  lack  of  predicted  dune  lowering  due  to  dune 
overwash  processes.  Figure  54  shows  the  measured  and  simulated  recession 
of  the  2.5  m  contour.  The  CSHORE  simulations  indicate  recession  of  the 
2.5  m  contour  at  less  than  5  m  at  profiles  37, 45,  56  and  103  whereas  the 
measurements  indicate  erosion  of  the  2.5  m  contour  at  10  to  20  m  at  the 
profiles  63  and  74.  Measured  and  computed  landward  storm  intrusion  are 
compared  in  Figure  55  which,  shows  that  CSHORE  predicts  a  landward 
storm  intrusion  that  compares  favorably  to  measurements  at  profiles  45,  74, 
and  103  but  does  not  simulate  the  landward  storm  intrusion  well  at  the 
other  profiles. 


Volume  Change  (m3/m) 
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Volume  Change  above  NGVD 
Ocean  City,  Maryland  (Nov-Jan  Storms) 


Figure  53.  Comparison  of  measured  and  simulated  volume  change  above  NGVD 
for  the  NJ  storm  series,  Ocean  City,  MD. 


Recession  of  +2.5  m  Contour 
Ocean  City,  Maryland  (Nov-Jan  Storms) 


Figure  54.  Comparison  of  measured  and  simulated  recession  of  the  2.5-m 
contour  for  the  NJ  storm  series,  Ocean  City,  MD. 
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Landward  Storm  Intrusion 
Ocean  City,  Maryland  (Nov-Jan  Storms) 


Figure  55.  Comparison  of  measured  and  simulated  storm  intrusion  for  the  NJ 
storm  series,  Ocean  City,  MD 

October- January  storms.  The  October-January  (OJ)  cases  involved 
the  simulation  of  the  entire  three  storm  sequence  for  four  profile  lines  that 
were  not  surveyed  in  November  1991  between  the  Halloween  and 
November  storms.  For  these  cases,  the  pre-storm  profiles  were  from  the 
June  1991  surveys  and  the  post-storm  profiles  were  obtained  from  the 
surveys  performed  after  the  January  storm.  Figure  56  shows  the  extent  to 
which  mass  was  conserved  between  the  pre-  and  post-storm  profiles. 
Although  mass  was  more  or  less  conserved  at  profiles  52,  81,  and  86, 
profile  92  gained  considerable  volume  between  the  surveys. 

The  CSHORE  simulation  results  for  the  OJ  storm  series  are  shown  in 
Figures  57  through  60.  These  figures  show  that  the  CSHORE  simulations 
provide  a  good  reproduction  of  the  measured  storm  induced  beach  profile 
change  in  the  absence  of  dune  overwash  processes.  The  measured  beach 
profile  response  for  profiles  52  and  92  are  reasonably  well  estimated 
whereas,  the  measured  dune  lowering  observed  at  profiles  81  and  86  is  not 
reproduced  by  the  CSHORE  predictions. 
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Mass  Conservation 

Ocean  City,  Maryland  (Oct -Jan  Storms) 


Figure  56.  Conservation  of  mass  between  beach  profile  surveys  for  the  OJ  storm  series, 

Ocean  City,  MD. 


Figure  57.  CSHORE  simulation  of  the  OJ  storm  series  for  profile  52,  Ocean  City,  MD. 
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Figure  58.  CSHORE  simulation  of  the  OJ  storm  series  for  profile  81  Ocean  City,  MD. 
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Figure  59.  CSHORE  simulation  of  the  OJ  storm  series  for  profile  86,  Ocean  City,  MD. 
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OJ92 

Oceancity  MD  profile  92  Oct  thru  Jan  (Calibrated  to  OJ81 B) 
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Figure  60.  CSHORE  simulation  of  the  OJ  storm  series  for  profile  92,  Ocean  City,  MD. 

The  quantitative  model  performance  measures  for  the  OJ  storm  series  cases 
are  illustrated  in  Figures  6i  through  63.  In  Figure  61  it  is  seen  that  the 
CSHORE  predicted  volume  change  above  NGVD  is  within  15  to  20  cu  m/m 
at  all  profiles  except  profile  81.  Measured  and  computed  recession  of  the 
2.5-m  contour  is  compared  in  Figure  62  where  excellent  agreement  is 
observed  at  profiles  52  and  92,  poorest  at  profile  81  and  agreement  within 
about  8  m  is  observed  at  profile  86.  The  predicted  and  measured  storm 
intrusion  distances  shown  in  Figure  63  are  within  about  10  m  at  profiles  52 
and  86,  less  than  20  m  at  profile  92,  and  just  over  20  m  at  profile  81.  It  is 
noted  that  the  measurements  show  a  considerable  alongshore  variation  in 
the  beach  profile  response  between  the  various  profiles  whereas  the  model 
results  are  more  consistent.  These  differences  can  be  at  least  partially 
explained  by  noting  that  the  environmental  forcing  (waves  and  water  levels) 
in  CSHORE  is  constant  between  the  profiles,  whereas  the  measurements 
reflect  non-uniform  beach  profile  response  that  arises  from  non-uniform 
forcing  due  to  known  irregularities  in  the  nearshore  bathymetry  offshore  of 
Ocean  City,  MD  (see  Stauble  et  al.  1993).  Specifically,  profiles  81  and  86,  are 
known  to  lie  in  an  erosional  hot  spot  area  that  is  believed  to  be  associated 
with  offshore  shore  oblique  shoals. 
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Volume  Change  above  NGVD 
Ocean  City,  Maryland  (Oct-Jan  Storms) 


Figure  61.  Comparison  of  measured  and  simulated  volume  change  above  NGVD  for  the 

OJ  storm  series,  Ocean  City,  MD. 


Recession  of  +2.5  m  Contour 
Ocean  City,  Maryland  (Oct-Jan  Storms) 


Figure  62.  Comparison  of  measured  and  simulated  recession  of  the  2.5-m 
contour  for  the  OJ  storm  series,  Ocean  City,  MD. 
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Landward  Storm  Intrusion 
Ocean  City,  Maryland  (Oct-Jan  Storms) 


Figure  63.  Comparison  of  measured  and  simulated  storm  intrusion  for  the  OJ 
storm  series,  Ocean  City,  MD. 

Manasquan  and  Point  Pleasant  Beach,  NJ 

Manasquan  and  Point  Pleasant  Beach  are  the  coastal  communities 
immediately  to  the  north  and  south  of  Manasquan  Inlet  on  the  Atlantic 
coast  of  New  Jersey.  The  beach  at  Manasquan  is  relatively  narrow  with 
short  groins  spaced  at  200  to  300-m  intervals.  Point  Pleasant  Beach  on 
the  other  hand  has  a  comparatively  wide  beach  with  no  littoral  structures 
other  than  the  south  jetty  at  Manasquan  Inlet.  These  contrasting  beach 
conditions  are  explained  by  the  net  northerly  longshore  sand  transport 
regime  and  the  inlet  stabilization  structures  at  Manasquan  Inlet.  Blockage 
of  the  net  longshore  transport  rate  by  the  jetties  at  Manasquan  Inlet 
produces  both  the  wide  beaches  in  Point  Pleasant  Beach  and  the 
comparatively  narrow  beaches  in  Manasquan. 

A  northeaster  impacted  these  two  beaches  between  28  and  30  March  1984 
causing  severe  erosion.  The  effects  of  this  storm  were  captured  by  profile 
surveys  conducted  by  the  U.S.  Army  Engineer  District  (USAED), 
Philadelphia,  as  part  of  the  Corps  of  Engineers  Monitoring  of  Completed 
Coastal  Projects  (MCCP)  program.  Subaerial  profile  surveys  at  the  two 
beaches  were  taken  approximately  two  days  before  the  storm  and  again 
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three  to  four  days  after  the  storm.  Subaqueous  profile  surveys  were  taken 
approximately  three  months  before  and  two  weeks  after  the  storm.  The 
subaerial  and  subaqueous  profiles  were  combined  to  obtain  an  approxi¬ 
mation  of  the  pre-  and  post-storm  profiles.  Wave  data  for  the  northeaster 
were  provided  by  a  waverider  buoy  located  at  a  depth  of  15.2  m  operated  as 
another  element  of  the  MCCP  program.  Water  level  data  were  obtained 
through  a  National  Oceanographic  Service  (NOS)  maintained  tide  gauge 
located  near  the  shoreward  side  of  Manasquan  Inlet.  The  time  histories  of 
the  wave  and  water  level  data  are  shown  in  Figure  64.  The  storm  duration 
was  approximately  2  days  with  a  peak  significant  wave  height  of  7  m  and  a 
peak  water  level  of  2  m  NGVD. 


Figure  64.  Wave  height,  wave  period,  and  water  elevation  time-histories  for  the  March  1984 
storm,  Manasquan/Point  Pleasant  Beach,  NJ. 
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Manasquan.  The  total  volume  change  per  unit  width  of  beach  based  on  the 
pre-  and  post-storm  beach  profiles  at  the  Manasquan  is  illustrated  in 
Figure  65.  The  beach  at  Manasquan  is  located  on  the  north  side  of 
Manasquan  Inlet.  The  beach  profiles  at  Manasquan  are  numbered  in  order 
with  distance  from  the  inlet  jetty  with  profile  Ml  being  the  closest  and  M9  the 
most  distant.  As  seen  in  Figure  65,  mass  conservation  between  the  pre-  and 
post-storm  beach  profiles  at  Manasquan  was  poor.  Five  of  the  nine  profiles 
were  found  to  have  either  gained  or  lost  more  than  100  cu  m/m.  These  large 
volumetric  changes  are  presumed  to  be  attributed  to  the  long  time  interval 
between  the  subaqueous  profile  surveys,  gradients  in  the  longshore  sand 
transport  due  to  the  inlet  structures,  or  possible  survey  error. 

CSHORE  simulation  results  for  the  March  1984  storm  at  Manasquan  are 
illustrated  in  Figures  66  to  74.  In  general,  the  CSHORE  predictions 
underestimate  the  measured  beach  profile  response  in  terms  of  total 
volumetric  erosion  particularly  erosion  high  on  the  beach  profile. 


Mass  Conservation 
Manasquan,  NJ 


Figure  65.  Conservation  of  mass  between  beach  profile  surveys  for  the  March  1984  storm, 

Manasquan,  NJ. 
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M3 

Profile:  M3  Manasquan 


s  68.  CSHORE  simulation  of  the  March  1984  storm  for  profile  M3,  Manasquan,  NJ 


M4 

Profile:  M4  Manasquan 


Figure  69.  CSHORE  simulation  of  the  March  1984  storm  for  profile  M4,  Manasquan,  NJ 


Elevation  [m]  NGVD  ^  |  Elevation  [m]  NGVD 


ERDC/CHL  TR-12-22 


ITT 


70.  CSHORE  simulation  of  the  March  1984  storm  for  profile  M5,  Manasquan,  NJ 


Figure  71.  CSHORE  simulation  of  the  March  1984  storm  for  profile  M6,  Manasquan,  NJ 
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M7 

Profile:  M7  Manasquan 


3  72.  CSHORE  simulation  of  the  March  1984  storm  for  profile  M7,  Manasquan,  NJ 


M8 

Profile:  M8  Manasquan 


Figure  73.  CSHORE  simulation  of  the  March  1984  storm  for  profile  M8,  Manasquan,  NJ 
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M9 

Profile:  M9  Manasquan 


Figure  74.  CSHORE  simulation  of  the  March  1984  storm  for  profile  M9,  Manasquan,  NJ. 

A  curious  characteristic  of  the  Manasquan  beach  profile  response  to  the 
March  1984  storm  is  the  lack  of  significant  erosion  between  the  elevations 
of  o  and  2  m  NVGD  with  an  associated  large  amount  of  erosion  above  the 
2  m  elevation.  The  specific  cause  of  this  unusual  beach  profile  response  is 
unknown  although  it  could  be  that  the  reduced  erosion  in  the  o  to  2  m 
elevation  range  is  the  result  of  the  rapid  development  of  a  post-storm 
recovery  berm  in  the  waning  stages  of  the  storm  and  prior  to  the  post¬ 
storm  profile  surveys. 

The  quantitative  model  performance  measures  for  the  March  1984  storm  at 
Manasquan  are  illustrated  in  Figures  75  through  77.  Figure  75  shows  that 
the  predicted  volume  loss  above  NGVD  is  consistently  under  predicted, 
except  at  profile  M7.  Predicted  recession  of  the  3.75  m  contour  is 
consistently  under  estimated  by  1  to  12  m  in  the  model  simulations,  as  seen 
in  Figure  76.  The  measurement  of  landward  storm  intrusion  illustrated  in 
Figure  77  shows  that  the  CSHORE  predictions  also  under  estimated  the 
storm  intrusion  distance  at  all  profiles  by  approximately  3  to  12  m. 
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Landward  Storm  Intrusion 
Manasquan,  NJ 


Profile 

Figure  77.  Comparison  of  measured  and  simulated  storm  intrusion  for  the  March  1984 

storm,  Manasquan,  NJ. 

Point  Pleasant  Beach.  Eight  profiles  were  surveyed  at  Point  Pleasant 
Beach  located  south  of  Manasquan  Inlet.  The  profiles  at  Point  Pleasant 
Beach  are  numbered  in  order  with  distance  from  the  inlet.  Profile  Pi  was 
closest  to  the  inlet  jetty  and  profile  P8  was  the  most  distant.  Figure  78 
shows  the  total  measured  volume  change  between  the  pre-  and  post-storm 
profile  surveys.  In  this  figure  it  is  seen  that  seven  of  the  eight  profiles 
experienced  a  net  loss  of  volume  between  the  pre-  and  post-storm  surveys. 
The  smaller  magnitude  of  the  net  volume  changes  at  the  Point  Pleasant 
Beach  profiles  compared  to  Manasquan  profiles  indicates  better  mass 
conservation.  The  Point  Pleasant  Beach  profiles  experienced  more  erosion 
and  recession  of  the  upper  profile  than  was  observed  on  the  Manasquan 
profiles.  Nearly  all  of  the  Point  Pleasant  Beach  post-storm  profiles 
contained  a  ridge  at  the  shoreline  which  is  likely  the  development  of  a 
post-storm  recovery  berm  formed  by  long  period  waves  in  the  waning 
stages  of  the  storm. 
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Mass  Conservation 
Point  Pleasant  Beach,  NJ 


Figure  78.  Conservation  of  mass  between  beach  profile  surveys  for  the  March  1984 
storm,  Point  Pleasant  Beach,  NJ. 

CSHORE  simulation  results  for  the  March  1984  storm  at  Point  Pleasant 
Beach  are  illustrated  in  Figures  79  through  86.  Overall,  the  CSHORE 
predictions  underestimate  the  measured  beach  profile  erosion  on  the 
beach  profile  above  the  2  m  NGVD  elevation.  Erosion  on  the  subaerial 
profile  between  the  o  and  2  m  NGVD  contours  is  over  estimated  in  the 
CSHORE  results.  However,  the  presence  of  a  post-storm  recovery  berm  in 
the  measured  profile  data  may  be  masking  the  storm  generated  erosion  in 
this  elevation  range. 

The  quantitative  model  performance  measures  for  the  March  1984  storm  at 
Point  Pleasant  Beach  are  illustrated  in  Figures  87  through  89.  Figure  87 
shows  that  the  predicted  volume  loss  above  NGVD  is  in  good  agreement 
with  the  measurements  (differences  are  less  than  10  cu  m/m)  for  five  of  the 
eight  profiles  although  the  measured  and  predicted  post-storm  profiles 
(Figures  79-86)  show  that  the  model  predictions  underestimated  erosion 
high  on  the  profile.  The  predictive  skill  implied  by  the  volume  loss  above 
NGVD  performance  measure  is  tempered  by  the  lack  of  predictive  skill 
illustrated  by  the  performance  measures  of  contour  recession  and  storm 
intrusion.  Figure  88  shows  that  recession  of  the  2.75  m  contour  is  dramati¬ 
cally  underestimated  in  the  CSHORE  predictions  for  six  of  the  eight  Point 
Pleasant  Beach  profiles.  Likewise,  the  extent  of  landward  storm  intrusion  is 
also  under  estimated  at  all  the  Point  Pleasant  Beach  profiles  except  Pi. 
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Profile:  PI  Point  Pleasant  (Default) 


;  79.  CSHORE  simulation  of  the  March  1984  storm  for  profile  PI,  Point  Pleasant 

Beach,  NJ. 


P2 

Profile:  P2  Point  Pleasant  (Default) 


Figure  80.  CSHORE  simulation  of  the  March  1984  storm  for  profile  P2,  Point  Pleasant 

Beach,  NJ. 
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;  81.  CSHORE  simulation  of  the  March  1984  storm  for  profile  P3,  Point  Pleasant 

Beach,  NJ. 
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Figure  82.  CSHORE  simulation  of  the  March  1984  storm  for  profile  P4,  Point  Pleasant 

Beach,  NJ. 


Elevation  [m|  NGVD  ^  Elevation  [m]  NGVD 


ERDC/CHL  TR-12-22 


P5 


0  100  200  300  400  500  600  700  800 

x|m] 


83.  CSHORE  simulation  of  the  March  1984  storm  for  profile  P5,  Point  Pleasant 

Beach,  NJ. 
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Figure  84.  CSHORE  simulation  of  the  March  1984  storm  for  profile  P6,  Point  Pleasant 

Beach.  NJ. 
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Figure  85.  CSHORE  simulation  of  the  March  1984  storm  for  profile  P7,  Point  Pleasant 

Beach,  NJ. 


Figure  86.  CSHORE  simulation  of  the  March  1984  storm  for  profile  P8,  Point  Pleasant 

Beach.  NJ. 
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Landward  Storm  Intrusion 
Point  Pleasant  Beach,  NJ 


Figure  89.  Comparison  of  measured  and  simulated  storm  intrusion  for  the 
March  1984  storm,  Point  Pleasant  Beach,  NJ. 

Dewey  Beach  and  Rehoboth  Beach,  NJ 

The  communities  of  Dewey  Beach  and  Rehoboth  Beach  on  the  east  coast  of 
Delaware  were  impacted  by  a  northeaster  storm  event  on  to  December 
1992.  The  beach  response  to  the  storm  was  captured  by  profile  surveys 
collected  by  the  USACE,  Philadelphia,  as  part  of  an  ongoing  feasibility 
study.  Pre-storm  profiles  were  surveyed  on  29  October  1992  and  post¬ 
storm  surveys  were  performed  after  passage  of  the  storm  on  18  December 
1992.  Wave  conditions  during  the  storm  were  recorded  by  a  wave  gauge 
located  at  a  depth  of  approximately  9  m  offshore  of  Dewey  Beach.  Tide 
information  was  available  from  a  nearby  tide  gauge  located  at  Lewes,  DE. 
The  time  histories  of  wave  and  water  level  information  for  the  10 
December  storm  event  are  plotted  in  Figure  90.  The  storm  had  a  peak 
significant  wave  height  of  4  m  and  a  peak  total  water  elevation  of  just  less 
than  2  m  NGVD  over  a  duration  of  approximately  four  days.  The  sediment 
grain  size  for  the  Dewey  Beach  and  Rehoboth  Beach  CSHORE  simulations 
was  specified  at  0.33  mm,  based  on  the  sediment  analyses  performed  as 
part  of  the  USAED,  Philadelphia  project  feasibility  study. 
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Figure  90.  Wave  height,  wave  period,  and  water  elevation  time-histories  for  the  December 
1992  storm,  Dewey  Beach/Rehoboth  Beach,  DE. 


Dewey  Beach.  Pre-  and  post-storm  profile  measurements  were  available 
for  seven  profile  lines  at  Dewey  Beach.  Figure  91  shows  the  storm  induced 
volume  change  between  the  pre-  and  post-storm  surveys  at  each  profile 
line.  The  figure  shows  that  a  net  volume  loss  was  experienced  at  all  of  the 
Dewey  Beach  profile  lines.  The  greatest  volume  loss  was  found  at  profile 
240  with  a  loss  of  approximately  120  cu  m/m  whereas  mass  was  nearly 
conserved  at  profile  140  where  the  volume  loss  was  less  than  10  cu  m/m. 
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Figure  91.  Conservation  of  mass  between  beach  profile  surveys  for  the  December  1992 

storm,  Dewey  Beach,  DE. 


The  CSHORE  simulations  of  storm  induced  beach  profile  change  at  Dewey 
Beach  provided  remarkably  good  agreement  with  the  measured  profile 
changes  as  illustrated  in  Figures  92  through  98.  Near  complete  erosion  of 
the  beach  berm  together  with  moderate  erosion  of  the  seaward  face  of  the 
dune  with  deposition  in  an  offshore  bar  feature  is  typical  of  measured  and 
predicted  beach  profile  response  to  the  December  1992  storm  at  Dewey 
Beach.  CSHORE  accurately  predicted  the  position  and  form  of  the  offshore 
bar  particularly  on  profiles  140,  210,  220,  and  240.  Erosion  of  the  upper 
beach  is  also  accurately  predicted  except  for  the  near  horizontal  swash 
platform  at  about  o  NGVD  evident  in  most  of  the  measured  profiles. 

The  quantitative  model  performance  measures  for  the  December  1992 
storm  at  Dewey  Beach  are  illustrated  in  Figures  99  through  101.  Figure  99 
shows  that  the  difference  between  the  measured  and  computed  volume 
change  above  o  NGVD  is  less  than  10  cu  m/m  for  two  of  the  seven  Dewey 
Beach  profiles  but  greater  than  or  equal  to  about  20  cu  m/m  for  the 
remainder.  The  predicted  recession  of  the  1.5-meter  contour  differs  from 
the  measured  by  less  than  10  m  in  all  Dewey  Beach  cases  (Figure  100).  The 
comparison  of  measured  and  simulated  storm  intrusion  distance  in 
Figure  101  also  shows  very  good  agreement  with  model  estimates  within 
5  m  of  the  measurements. 
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92.  CSHORE  simulation  of  the  December  1992  storm  for  profile  D100,  Dewey 

Beach,  DE. 
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Figure  93.  CSHORE  simulation  of  the  December  1992  storm  for  profile  D115,  Dewey 

Beach,  DE. 
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i  94.  CSHORE  simulation  of  the  December  1992  storm  for  profile  D125,  Dewey 

Beach,  DE. 
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Figure  95.  CSHORE  simulation  of  the  December  1992  storm  for  profile  D140,  Dewey 

Beach,  DE. 
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3  96.  CSHORE  simulation  of  the  December  1992  storm  for  profile  D210,  Dewey 

Beach,  DE. 
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Figure  97.  CSHORE  simulation  of  the  December  1992  storm  for  profile  D220,  Dewey 

Beach,  DE. 
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Figure  98.  CSHORE  simulation  of  the  December  1992  storm  for  profile  D240,  Dewey 

Beach,  DE. 


Figure  99.  Comparison  of  measured  and  simulated  volume  change  above  NGVD  for  the 
December  1992  storm,  Dewey  Beach,  DE. 
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re  100.  Comparison  of  measured  and  simulated  recession  of  the  1.5-m 
contour  for  the  December  1992  storm,  Dewey  Beach,  NJ. 
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Figure  101.  Comparison  of  measured  and  simulated  storm  intrusion  for  the 
December  1992  storm,  Dewey  Beach,  NJ. 


Dewey  Beach,  DE 


ERDC/CHL  TR-12-22 


119 


Rehoboth  Beach.  Rehoboth  Beach  lies  approximately  1.5  miles  north  of 
Dewey  Beach  along  the  Delaware  coast.  Pre-  and  post-storm  profile 
measurements  at  Rehoboth  Beach  were  available  at  eight  profile  lines  to 
model  the  beach  response  to  the  December  1992  storm.  Figure  102  shows 
the  storm  induced  volume  change  between  the  pre-  and  post-storm 
surveys  for  each  of  the  Rehoboth  Beach  profiles.  Mass  was  reasonably 
conserved  at  four  of  the  Rehoboth  Beach  profiles  (117,  200,  215,  and  224) 
whereas  net  volumetric  losses  of  40  to  60  cu  m/m  were  experienced  at  the 
other  four  profiles  (122, 131, 138,  and  236). 


Mass  Conservation 
Rehoboth  Beach,  DE 


Profile 

Figure  102.  Conservation  of  mass  between  beach  profile  surveys  for  the  December 
1992  storm,  Rehoboth  Beach,  DE. 

The  CSHORE  simulation  results  for  the  December  1992  storm  at  Rehoboth 
Beach,  DE  are  illustrated  in  Figures  103  through  110.  These  figures  show  the 
presence  of  an  approximately  50  m  wide  berm  backed  by  a  dune  or  seawall 
in  the  pre-storm  profiles.  The  CSHORE  simulations  predicted  complete 
erosion  of  the  berm  feature  together  with  slight  erosion  low  on  the  seaward 
face  of  the  dune  in  good  agreement  with  the  measured  post-storm  profiles. 
Profiles  117, 122,  and  200  of  the  Rehoboth  Beach  data  set  included  a  seawall 
landward  of  the  beach  berm.  Although  CSHORE  does  not  presently  include 
a  seawall  boundary  condition  the  predicted  beach  profile  response  is 
excellent  seaward  of  the  seawall  position. 
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Figure  103.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R117, 

Rehoboth  Beach,  DE. 


Figure  104.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R122, 

Rehoboth  Beach,  DE. 
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Figure  105.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R131, 

Rehoboth  Beach,  DE. 


Figure  106.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R138, 

Rehoboth  Beach,  DE. 
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107.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R200, 
Rehoboth  Beach,  DE. 


Figure  108.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R214, 

Rehoboth  Beach,  DE. 
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Figure  109.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R224, 

Rehoboth  Beach,  DE. 


Figure  110.  CSHORE  simulation  of  the  December  1992  storm  for  profile  R236, 

Rehoboth  Beach,  DE. 
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The  quantitative  model  performance  measures  for  the  December  1992 
storm  at  Rehoboth  Beach  are  illustrated  in  Figures  111  through  113. 

Figure  111  shows  that  overall  the  CSHORE  predictions  of  volumetric  change 
above  NGVD  compare  well  with  measured  values  with  the  exception  of 
profile  236.  In  five  of  the  profiles  the  volumetric  change  estimate  was  within 
about  10  cu  m/m,  profile  R117  was  within  15  cu  m/m,  and  profile  R236  was 
under  estimated  by  approximately  40  cu  m/m.  Some  of  the  overestimation 
of  volume  change  above  NGVD  could  be  attributed  to  the  presence  of  a  post¬ 
storm  recovery  berm  in  the  measured  post-storm  profiles.  Figure  112  shows 
the  measured  and  calculated  recession  of  2  m  contour  where  it  is  seen  that 
the  difference  between  the  predicted  and  measured  recession  is  less  than 
5  m  for  all  but  one  of  the  profiles  and  approximately  9  m  at  the  other  profile. 
Figure  113  shows  good  agreement  between  the  measured  and  calculated 
landward  storm  intrusion  distance.  The  measured  and  computed  storm 
intrusion  distances  for  the  Rehoboth  Beach  simulations  differ  by  less  than 
10  m  for  five  of  the  eight  profiles.  The  greatest  difference  is  approximately 
18  m  at  profile  236. 
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Figure  111.  Comparison  of  measured  and  simulated  volume  change  above  NGVD  for  the 
December  1992  storm,  Dewey  Beach,  DE. 
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Figure  112.  Comparison  of  measured  and  simulated  recession  of  the  1.5-m 
contour  for  the  December  1992  storm,  Dewey  Beach,  NJ. 
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Figure  113.  Comparison  of  measured  and  simulated  storm  intrusion  for  the 
December  1992  storm,  Dewey  Beach,  NJ. 
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Debidue  Beach  and  Myrtle  Beach,  SC 

Hurricane  Hugo  made  landfall  near  Charleston,  SC  on  22  September  1989. 
The  beaches  along  the  South  Carolina  coast  north  of  Charleston  were 
exposed  to  strong  onshore  winds  and  high  waves  and  water  levels.  The 
extreme  erosional  effects  of  Hurricane  Hugo  were  capture  by  beach 
surveys  performed  by  Coastal  Science  and  Engineering  Inc.  and  the  U.S. 
Army  Corps  of  Engineers,  Waterways  Experiment  Station,  Coastal 
Engineering  Research  Center  (Stauble  et  al.  1990).  Beach  profile  data 
collected  at  Debidue  Beach  and  Myrtle  Beach  in  June  1989  provided  the 
pre-storm  profiles  for  modeling  the  beach  response  to  Hurricane  Hugo.  A 
single  post-storm  profile  collected  within  a  week  of  the  passing  of  the 
hurricane  was  available  for  simulating  the  impact  of  the  hurricane  at 
Debidue  Beach  and  Myrtle  Beach.  A  median  sediment  grain  size  of 
0.20  mm  was  specified  in  the  model  simulations  as  a  representative  value 
for  both  beaches. 

Water  surface  elevations  during  Hurricane  Hugo  were  recorded  by  a  tide 
gauge  located  in  Winyah  Bay,  approximately  10  miles  south  of  Debidue 
Beach  and  40  miles  south  of  Myrtle  Beach.  These  tide  gauge  data 
represented  the  best  available  estimates  of  the  total  water  level  variations 
produced  by  the  hurricane  at  both  beaches.  It  is  believed  that  the  actual 
storm  surge  at  the  study  sites  was  probably  somewhat  lower  than  the 
measurements  since  the  beaches  were  located  farther  north  away  from  the 
center  of  the  storm.  Hindcast  wave  information  was  used  as  input  for  the 
CSHORE  simulation  because  no  local  measurements  were  available.  To 
account  for  spatial  variation  in  wave  conditions,  wave  information  from 
different  hindcast  stations  were  used  for  the  two  beaches.  The  time 
histories  of  the  hindcast  wave  height  and  period  and  measured  water  level 
used  in  the  modeling  are  shown  in  Figures  114  and  115.  As  indicated  in  the 
figures  the  peak  hindcast  significant  wave  height  at  Debidue  Beach  was 
approximately  6  m  whereas  at  Myrtle  Beach  the  peak  hindcast  significant 
wave  height  was  approximately  4  m.  The  differences  in  the  hindcast  wave 
conditions  at  the  two  sites  are  primarily  related  to  the  distance  from  the 
storm  center  to  the  subject  beach.  The  maximum  water  elevation 
measured  at  the  Winyah  Bay  tide  gauge  was  just  less  than  2.5  m  NGVD. 
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Figure  114.  Wave  height,  wave  period  and  water  elevation  time-histories  for  Hurricane 

Hugo,  Debidue  Beach,  SC. 
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Figure  115.  Wave  height,  wave  period  and  water  elevation  time-histories  for  Hurricane 

Hugo,  Myrtle  Beach,  SC. 


CSHORE  simulation  results  for  Hurricane  Hugo  at  Debidue  Beach  and 
Myrtle  Beach  are  presented  in  Figures  116  and  117.  Figure  116  shows  that 
CSHORE  significantly  under  predicted  the  near  complete  erosion  of  the 
berm  and  rather  wide  foredune  indicated  in  the  post-storm  measured 
profile.  The  measured  post-storm  profile  also  indicates  considerably  more 
nearshore  deposition  between  -4  m  and  +1  m  NGVD  than  predicted  by 
CSHORE.  Figure  117  shows  that  the  measured  upper  beach  profile  erosion 
due  to  Hurricane  Hugo  at  Myrtle  Beach  is  significantly  under  predicted  in 
the  CSHORE  simulation. 


Elevation  [m]  NGVD  Elevation  (m]  NGVD 
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Debidue  Beach-Hurricane  Hugo  (Default) 


Figure  116.  CSHORE  simulation  of  Hurricane  Hugo,  Debidue  Beach,  SC. 


MYR 

Myrtle  Beach-Hurricane  Hugo  (Default) 


Figure  117.  CSHORE  simulation  of  Hurricane  Hugo,  Myrtle  Beach,  SC. 
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The  quantitative  model  performance  measures  for  Hurricane  Hugo  and 
Debidue  Beach  and  Myrtle  Beach  are  presented  in  Figures  118  through 
120.  Figure  118  shows  that  the  CSHORE  calculated  erosion  volume  above 
NGVD  was  approximately  50  cu  m/m  and  33  cu  m/m  less  than  indicated 
by  the  pre-  and  post-storm  profile  measurements.  Recession  of  the  3.5  m 
NGVD  contour  presented  in  Figure  119  was  under  estimated  by 
approximately  20  m  at  Debidue  Beach  and  7  m  at  Myrtle  Beach. 

The  landward  storm  intrusion  distance  was  also  under  estimated  by  8  m 
and  5  m  as  indicated  in  Figure  120. 

Southern  California  Field  Cases 

On  January  18, 1988  a  significant  coastal  storm  impacted  the  beaches  of 
southern  California.  Wave  and  water  level  information  for  this  storm  was 
obtained  from  measurements  at  Wave  Gauge  051  of  the  Coastal  Data 
Information  Program  (CDIP,  http:/ /cdip.ucsd.edu/)  operated  and 
maintained  by  Scripps  Institute  of  Oceanography,  University  of  California, 
San  Diego.  This  storm  event  was  the  largest  recorded  wave  event  at  Wave 
Gauge  051  in  the  interval  September  1, 1983  through  June  13, 1988.  A  peak 
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Figure  118.  Comparison  of  measured  and  simulated  volume  change  above  NGVD  for 
Hurricane  Hugo,  Debidue  Beach  and  Myrtle  Beach,  SC. 
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Recession  of  +3.5  m  Contour 
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e  119.  Comparison  of  measured  and  simulated  recession  of  the  3.5-m 
contour  for  Hurricane  Hugo,  Debidue  Beach  and  Myrtle  Beach,  SC. 
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Figure  120.  Comparison  of  measured  and  simulated  storm  intrusion  for 
Hurricane  Hugo,  Debidue  Beach  and  Myrtle  Beach,  SC. 
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significant  wave  height  of  4.52  m  was  recorded  on  January  18, 1988,  which 
is  nearly  five  times  larger  than  the  mean  recorded  significant  wave  height  of 
0.92  m  for  the  previous  five  years.  The  time  histories  of  the  measured 
significant  wave  height,  peak  wave  period,  and  total  water  elevation  are 
plotted  in  Figure  121.  Here  it  is  seen  that  although  the  storm  generated  very 
large  wave  conditions  there  is  no  statistically  significant  rise  in  the  mean 
water  level  relative  to  the  tidal  variation  of  about  2  m.  Pre-  and  post-storm 
beach  profile  surveys  are  available  at  Oceanside  and  Del  Mar,  CA  for 
modeling  the  beach  profile  response  to  the  January  1988  storm.  The  pre¬ 
storm  profile  survey  was  obtained  in  September  1987  and  the  post-storm 
profile  survey  was  conducted  soon  after  passage  of  the  storm  in  late  January 
1988.  Mean  sediment  grain  sizes  of  0.17  mm  and  0.18  mm  were  specified 
for  the  Oceanside  and  Del  Mar  profiles,  respectively  based  on  sediment 
sampling. 

CSHORE  simulation  results  for  the  January  1988  storm  at  Oceanside  and 
Del  Mar,  CA  are  illustrated  in  Figures  122  and  123.  As  seen  in  Figurei20 
both  the  post-storm  measured  profile  and  the  CSHORE  prediction  indicate 
that  the  storm  produced  complete  erosion  of  the  berm  at  Oceanside,  CA. 

The  CSHORE  simulation  predicts  development  of  an  offshore  bar  feature  in 
agreement  with  the  measured  profile  but  under  predicts  the  development  of 
a  trough  landward  of  the  bar  and  erosion  of  the  forshore  slope  between  the  - 
4  m  contour  and  MLLW  shoreline.  At  Del  Mar  (Figure  123)  storm  induced 
erosion  above  the  1  m  contour  is  well  predicted  as  is  the  development  and 
position  of  the  offshore  bar  feature.  As  in  the  Oceanside  case  the  model  does 
not  predict  development  of  a  trough  feature  landward  of  the  bar. 

The  quantitative  model  performance  measures  for  the  Oceanside  and 
Del  Mar,  CA  field  cases  are  presented  in  Figures  124  through  126.  In 
Figure  124  it  is  seen  that  the  measured  volumetric  erosion  above  MLLW 
varies  from  approximately  95  cu  m/m  at  Oceanside  to  about  too  cu  m/m  at 
Del  Mar.  The  CSHORE  predicted  volumetric  erosion  above  MLLW  is  within 
1  cu  m/m  of  measurements  at  Del  Mar  and  approximately  13  cu  m/m  at 
Oceanside,  and  is  viewed  as  good  agreement.  Measured  and  computed 
recession  of  the  +3.0  m  contour  is  compared  in  Figure  125  where  it  is  seen 
that  the  model  predictions  are  within  about  1  or  2  m  of  the  measurements. 
The  difference  between  the  measured  and  computed  landward  storm 
intrusion  distance  is  seen  in  Figure  126  to  be  approximately  2  m  at 
Oceanside  and  5  m  at  Del  Mar.  Overall,  the  Oceanside  and  Del  Mar 
California  beach  profile  response  to  the  January  1988  storm  event  is  well 
predicted  by  CSHORE  using  the  west  coast  default  parameter  values. 
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Figure  122.  CSHORE  simulation  of  January  1988  storm,  Oceanside,  CA. 


DM580 
Del  Mar,  CA 

I  I  I  I  I  I 

1  Initial  Profile 
'  Measured  Final 

CSHORE,  Defaults,  RMS  error  =  0.27  m| 


effb=  0.0100 
eff,  =  0.0200 

a  =  0.5000 

b  =  0.0020 


100  200  300  400  500  600  700  800  900 

x  [m] 


Figure  123.  CSHORE  simulation  of  January  1988  storm,  Del  Mar,  CA. 
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Figure  124.  Comparison  of  measured  and  simulated  volume  change  above  MLLW  for  the 
January  1988  storm,  Oceanside  and  Del  Mar,  CA. 
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Figure  125.  Comparison  of  measured  and  simulated  recession  of  the  3.0-m 
contour  for  the  January  1988  storm,  Oceanside  and  Del  Mar,  CA. 
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Landward  Storm  Intrusion 
Oceanside  &  Del  Mar,  CA 


Figure  126.  Comparison  of  measured  and  simulated  storm  intrusion  for  the 
January  1988  storm,  Oceanside  and  Del  Mar,  CA. 
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11  Conclusions 

The  cross-shore  model  CSHORE  is  presented,  including  formulation 
documentation  and  validation  test  cases.  The  combined  wave  current 
model  CSHORE  based  on  the  time-averaged  continuity,  cross-shore 
momentum,  longshore  momentum,  wave  action  and  roller  energy 
equations  predicts  the  cross-shore  variations  of  the  mean  and  standard 
deviation  of  the  free  surface  elevation  and  depth-averaged  cross-shore  and 
longshore  velocities  under  normally  or  obliquely  incident  irregular 
breaking  waves.  The  sediment  transport  formulas  for  the  cross-shore  and 
longshore  transport  rates  of  suspended  sediment  and  bedload  on  sand 
beaches  are  relatively  simple  and  require  the  hydrodynamic  input 
variables  which  can  be  predicted  efficiently  and  fairly  accurately  using 
existing  wave  and  current  models.  The  CSHORE  model  has  been 
compared  with  a  number  of  small-scale  and  large-scale  laboratory  data 
and  is  ready  for  practical  applications.  CSHORE  has  been  extended  to  the 
intermittently  wet  and  dry  zone  for  the  prediction  of  wave  overwash,  levee 
erosion  and  deformation  of  a  low-crested  stone  structure. 

The  beach  erosion  predictive  capability  of  the  CSHORE  numerical 
simulation  model  was  evaluated  in  field  case  studies  using  data  describing 
severe  beach  erosion  at  nine  different  sites  for  eight  storms.  Seven  of  the 
field  sites  were  located  on  the  mid- Atlantic  east  coast  (South  Carolina  to 
New  Jersey)  and  two  were  located  on  the  southern  California  coast.  Model 
performance  focused  on  the  capability  of  the  model  to  predict  erosion  of  the 
beach  berm  and  dune.  Model  performance  was  quantified  by  comparison  of 
measured  and  calculated  volume  change,  recession  of  select  profile  con¬ 
tours,  and  the  landward  progression  of  damaging  erosion  (0.3  m  vertical 
erosion).  CSHORE  was  applied  to  the  field  data  sets  using  default  parameter 
values  estimated  based  on  site-specific  calibration  of  the  model  at  three  of 
the  field  sites.  Overall  the  model  provided  satisfactory  estimates  of  storm 
induced  beach  erosion  for  various  pre-storm  beach  profile  shapes  and 
conditions.  These  field  applications  of  CSHORE  demonstrated  model 
effectiveness  in  simulating  beach  erosion  but  also  identified  areas  requiring 
research  and  model  refinement.  Dune  lowering  and  landward  translation  of 
the  dune  crest  due  to  overwash  processes  during  the  January  4, 1991  storm 
at  Ocean  City,  MD  was  under  predicted  by  CSHORE.  The  magnitude  of 
upper  beach  erosion  was  also  under  predicted  for  Hurricane  Hugo  at 
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Debidue  Beach  and  Myrtle  Beach,  SC  as  well  as  at  six  of  the  eight  profiles  at 
Point  Pleasant  Beach,  NJ  due  to  a  northeaster  storm  event  in  March  1984. 
Improvements  are  needed  in  the  area  of  sediment  transport  and 
morphology  evolution  in  the  intermittently  wet  zone. 
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